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ABSTRACT 

We present an analysis of stellar population gradients in 4,546 Early-Type Galax- 
ies (ETGs) with photometry in grizYHJK along with optical spectroscopy. ETGs 
were selected as bulge-dominated systems, displaying passive spectra within the SDSS 
fibers. A new approach is described which utilizes color information to constrain age 
and metallicity gradients. Defining an effective color gradient, V^^, which incorporates 
all of the available color indices, we investigate how V^, varies with galaxy mass prox- 
ies, i.e. velocity dispersion, stellar (A/*) and dynamical (M^yn) masses, as well as age, 
metallicity, and [a/Fe]. ETGs with Mdyn larger than 8.5 x 10^° M© have increasing age 
gradients and decreasing metallicity gradients wrt mass, metallicity, and enhancement. 
We find that velocity dispersion and [a/Fe] are the main drivers of these correlations. 
ETGs with 2.5 x 10'" Mq < Mdyn < 8.5 x 10'° Mq, show no correlation of age, metal- 
licity, and color gradients wrt mass, although color gradients still correlate with stellar 
population parameters, and these correlations are independent of each other. In both 
mass regimes, the striking anti-correlation between color gradient and a-enhancement 
is significant at ~ 5 cr, and results from the fact that metallicity gradient decreases with 
[a/Fe]. This anti-correlation may reflect the fact that star formation and metallicity 
enrichment are regulated by the interplay between the energy input from supemovae, 
and the temperature and pressure of the hot X-ray gas in ETGs. For all mass ranges, 
positive age gradients are associated with old galaxies (> 5 - 7 Gyr). For galaxies 
younger than ~ 5 Gyr, mostly at low-mass, the age gradient tends to be anti-correlated 
with the Age parameter, with more positive gradients at younger ages. 
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Subject headings: Galaxies: formation - Galaxies: evolution- Galaxies: fundamental 
parameters 



Introduction 



The study of the stellar population in galaxies presents some unique challenges. First, the 
method used to characterize a stellar population is crucial as it provides constraints on theories 
of galaxy formation (e.g. De Lucia et al. 2006). Second, age and metallicity dominated indica- 
tors should minimize the age-metallicity degeneracy described in the 3/2 Worthey law (Worthey, 
Trager & Faber 1995). Finally, a galaxy's environment may be the dominant influence on its star 
formation history, further reducing the adequacy of the Single Stellar Population (SSP) assump- 
tion. Monolithic collapse models predict early-type systems with a central region that is more 
metal rich, slightly younger, and less a/F^ enhanced than the outer parts. Mergers may wash out 
any possible gradient present in the progenitors, although dissipation can trigger star formation in 
the cores of galaxies, creating age gradients. I t is importan t to emphasize that we still lack firm 
quantitative estimates for many of these effects (IWhitdl 19801) . Considering the impact of such pre- 
dictions, we need large and well characterized samples to tackle this issue. The nearby universe 
(z < 0.1) provides a starting point, as we now have high quality data for large numbers of galaxies. 

The study of color gradients is motivated by the possibility of discriminating between different 
models of galaxy formation, and points to the real need for systematic and accurate measurements 
for a large and well controlled sample of ETGs. These systems become gradually bluer in the 
outskirts (e.g. Peletier et al. 1990; Goudfrooij et al. 1994; Michard 1999; Wu et al. 2005; Cantiello 
et al. 2005; Suh et al. 2010). Although the age-metallicity degeneracy makes the interpretation of 
the color gradient troublesome, metallicity gradient seems to be main cause (e.g. Saglia et al. 2000; 
Tamura et al.ll2000l ; iTamura & Ohtall200Q ; La Barbera et al. 2003). Small, positive age gradients 
are also consistent with obs ervations (e.g. ISaglia et al.ll2000t) , and have been detec ted in the low- 
redshift population of ETGs ( La Barbera & de Carvalholl2009 ; Clemens et al.ll2009 ). Positive color 
gradients in high-redshift ETGs have been observed by lFerreras et al.l(|2005r) . who found that about 
one-third of field ETGs at z ~ 0.7 have blue cores, in contrast to only 10% at lower redshift (see 



also iMenanteau et al.ll200lLl2004l ; lFerreras et al.ll2009l) . The possibility that color gradients could 
be due merely to dust in a well mixed stellar population has been considered (e.g. Goudfrooij 
& de Jong 1995) but many ETGs exhibit radial gradients of metal absorption features, which 
cannot be attributed to dust (e.g. Gonzalez 1993; Mehlert et al. 2000). Furthermore, many ETGs 
are undetected by IRAS at 60 fi and 100 ju. Cosmological simulations do not yet have sufficient 
spatial resolution to provide reliable color gradients. However, chemo-dynamical simulations of 
ETGs are starting to yield meaningful results for metallicity gradients, reproducing some observed 
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population gradients (see Kobayashi 2004). 



The last decade has witnessed the appearance of important observational studies us ing photo 
metri c and spectroscopic dat a to constrain age, metallic i ty and elemental ab undances ( I J0rgensen 



19991 : IPoggianti et al. 



200 lal : iTerlevich & ForbesI l2002l : iTrager et all 120001) . They conclude that 



age and metallicity are anti-correlated for galaxies of a given luminosity. Furthermore, Poggianti 
et al. (2001b) find that in high density regimes a large fraction of systems shows no signs of star 
formation in their center s. Another impo r tant fi nding relating star formation and environment 
comes from the work of IKuntschner et al.l (12002) who analyzed a sample of early-type systems 
in the low density regime (field) showing that these galaxies are younger and more metal-rich 
than their counterparts in clusters. Most of these studies used samples of at most a few hundred 
galaxies. Since these works are the nearby counterparts of t he studies of the evolution of e arly- 
type galaxies (ETGs) as a fun ction of redshift in clusters (e.g. Ivan Dokkum & Stanfordll2003b and 
in low-density environments (|Treu et al.lll999h . it is necessary to establish large and well-defined 
samples covering a broad range of environments as we do here. 

This is the fourth paper of a series studying the global properties of ETGs in the local Universe 
(z < 0.1). We focus on stellar population gradients and how they depend on galaxy properties. We 
use opti cal data (^riz photometry and sp ectra) from the Sloan Digital Sky Survey Data Release 6 
(SDSS, lAdelman-McCarthy et al.ll2008[) and ne ar-infrared photometry from the UKIRT Infrared 
Deep Sky Survey Data Release Four (UKIDSS. Lawrence et al.ll2007|) . The selection of ETGs is 
described in La Barbera et al. (2010a; hereafter Paper I), and we refer the reader to that paper for 
all the details on sample selection and the procedures used to derive the galaxy parameters. The 
results presented in this paper apply not only to ellipticals, but also to lenticular galaxies as well as 
those dominated by a red central bulge (but see Sec. 3.3 and Sec. 8.3). 

This paper is organized as follows: Sec. 2 briefly describes the sample used while in Sec. 
3 we present how stellar and dynamical masses were computed and stellar population properties 
measured. Measurement of color gradients are described in Sec. 4, while their distributions are 
discussed in Sec. 5. Sec 6 details how we derive stellar population gradients from color gradients. 
Sec. 7 presents our approach to constrain the variation of stellar population parameters, and dis- 
cusses the systematics involved. Sec. 8 illustrates the dependence of the efl^ective color gradient on 
galaxy parameters and Sec. 9 examines the age versus metallicity gradients in ETGs. Finally, in 
Sec. 10 we summarize our most important results. Throughout the paper, we adopt a cosmology 
with Ho = 75 km s"^ Mpc^^ = 0.3, and Qa = 0.7. 
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2. The sample 



The sample of ETGs is selected from SDSS-DR6, following the procedure described in La 
Barbera et al. (2008b), La Barbera & de Carvalho (2009), and Paper I. We select galaxies in the 
redshift range of 0.05 to 0.095, with "^Mr< - 20, where °^M,. is the k-corrected SDSS Petrosian 
magnitude in r-band . The k-correction is determined using the software kcorrect (version AAA; 
Blanton et al.ll2003a ), through a restf rame r-band filter blue-shifted by a factor (1 -i- zo)- We adopt 
Zo = 0.1 as in (e. g.) iHogg et al.ll2004l The lower redshift limit of the sample is chosen to minimize 
the aperture bias (|G6mez et al. 2003 ), while the uppe r redshift limit guarantees not only a high level 
of completeness (according to lSorrentino et al. but also allows us to define a volume-limited 

sample of bright early-type systems. At the upper redshift limit of z = 0.095, the magnitude cut 
of -20 also corresponds approximately to the magnitude limit to whi ch the SDSS spectros copy 
is complete (i.e. a Petrosian magnitude of ~ 17.8). Following iBemardi et al.l (|2003h . we 
define ETGs using the SDSS spectroscopic parameter eClass, which indicates the spectral type 
of a galaxy on the basis of a principal component analysis, and the SDSS photometric parameter 
fracDev,., which measures the fraction of the galaxy light that is better fitted by a de Vaucouleurs 
(rather than an exponential) law. ETGs are those systems with eClass < and fracDeVr > 0.8. We 
select only galaxies with central velocity dispersion, ctq, available from SDSS-DR6, in the range 
of 70 and 420 km s^\ and with no spectroscopic warning on (i.e. zWarning attribute set to zero). 
These constraints are chosen to select only reliable velocity dispersion measurements from SDSS. 
Applying the above criteria leads to a sample of 39, 993 ETGs. 

Out of the 39, 993 ETGs defining the SPIDER sample, 5, 080 also have photometry in the 
YJHK wavebands from the UKIDSS-Large Area Survey (see Paper I). All galaxies have two 
measures of the central velocity dispersion, one from SDSS-D R6 and an alternate estirn ate ob- 
tained by fitting SDSS spectra with the STARLIGHT software (icid Femandes et al.ll2005 !). In all 
wavebands, the structural parameters including effective radius, r^, the mean surface brightness 
within tha t radius, < u >e, and th e Sersic index, n, have been homogeneously measured using 
2DPH0T (|La Barbera et al.ll2008ar) . Total magnitudes are computed from the effective parameters 
in each filter. 



For this work, we select ETGs from the SPIDER survey as follows. First, we exclude galaxies 
which are poorly fitted by a two-dimensional Sersic model. This is done by selecting only objects 
with < 2 in all wavebands, where is the rms of residuals between the galaxy image in 
a given band and the corresponding best-fitting two-dimensional Sersic model. The threshold 
= excludes those objects which are more than 5 <t above the peak of the;^^^ distribution in all 
passbands. As shown in Paper I, the distributions of;^^ in grizYJHK are peaked around one, with 
a width of ~ 0.2. Second, since we estimate the internal color gradients from and n, we exclude 
galaxies with large uncertainties on these parameters. As shown in Paper I, the errors on n are 
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nearly independent of S /N ratio, and are always less than ~ 0.2 dex. On the other hand, the errors 
on Tg strongly depend on the S IN. To remove galaxies with large uncertainties without overly 
reducing the sample size, we select only those ETGs with a log rg uncertainty less than 0.5 dex in 
all wavebands. These restrictions yield a subsample of 37, 068 galaxies with parameters in griz. 
Of these, 4, 546 objects also have photometry in YJHK. Hereafter, we refer to these as the optical 
and optical+NIR samples of ETGs, respectively. 



3. Galaxy parameters 

We analyze the dependence of internal color gradients in ETGs on galaxy parameters, i.e. 
structural parameters, optical and NIR magnitudes, stellar and dynamical masses, and stellar pop- 
ulation properties (age, metallicity, and a-enhancement). The derivation of magnitudes and struc- 
tural parameters is described in Paper I. Below, we describe the estimate of stellar and dynamical 
masses (Sec. l3.1l and l3.2l) . as well as that of the stellar population parameters (Sec. l3.3|) . 



3.1. Stellar masses from SED fitting 



Galaxy stella r masses can be estirnated b y fitting their observed SEDs with stellar population 
synthesis models. Brinchmann & Ellis ( 2000 f) demonstrated that such techniques, combining op- 
tical and near-IR fluxes with spectroscopic red shifts, could be u sed t o derive masses w ith a factor 



of two uncertainty to z ~ 1. Local studies by IbcU et all J2OO3I) and 



the stellar mass function using similar methods. Sali m et al. 



Co 



e et al.l (|2001 l) examined 



(l2005Ll2007h have shown that galaxy 



properties derived from SED fitting with broadband colors are consistent with those derived from 
spectroscopic data. Thus, we use SED fitting to derive galaxy stellar masses. 



We employ the Le Phare code (S. Amouts & O . Ilbert) used by the COSMOS survey (|Ilbert et al 



20091) to fit a suite of iBruzual & Chariot! (120031 hereafter BC03) model SEDs to our observed 
8-band photometry and the known spectroscopic redshift. LePhare performs a minimization 
between the observed fluxes (and their associated errors) and those synthesized from the models. 
B ecause early-types are unli kely to have significant ongoing star formation, we use a simple set 
of iBruzual & Chariot! (l2003b SEDs, with 6 difi'erent star-formation e-folding times r, 5 internal 
reddenings E{B - V), and four metallicities Z/Zp, for a total of 1 20 po ssible models. The values of 
these p arameters are given in Tab.[TJ We assume a lCalzetti et al. hood) extinction law andlChabrien 
(120031) initial mass function. The reddening E(B-V) is limited to 0.5 magnitudes to avoid incorrect 
fitting of observationally red galaxies as highly reddened blue galaxies; such objects should be 
absent from our sample. Because each BC03 template is normalized to one solar mass, we use the 
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observed A'-band flux to determine the normalization of the best-fit template which then gives the 
stellar mass. We use ^-band for the mass scaling since it is only we akly affected by dust extinc 



tion and is quite in sensitive to the presence of young, luminous stars (iKauffmann & Charlotlll998 



BundvetalJl2005 ). 



Galaxy fluxes are estimated in an adaptive aperture of radius 3 x vk^i in all wavebands, where 
r/f , is the i-band Kron radius (see Paper I). We use these aperture fluxes as they have a better 
accuracy than those obtained from (2DPH0T) total magnitudes. In order to derive total stellar 
masses, M^, we apply an aperture correction. We calculate the fraction of total to 3 x rjf , aperture 
flux in K-band, and then multiply the estimated stellar mass by this fraction. This procedure 
assumes that stellar mass in a galaxy is distributed in the same way as the K-band light. 

Rather than use an individual best-fitting SED (based on the x^)^ we use a mass computed 
from the model corresponding to the median e"^^^^. This reduces stochastic mass errors due to 
individual models which happen to have very low;^^^ values, since there remains some degeneracy 
in the SED fits even with 8 filters. A small number of galaxies (1465 of 39993 with optical data 
only, and 286 of 4546 with NIR data) are not well fit by any of the SEDs, resulting in = 

and thus no stellar mass estimate. These objects, comprising 3.6% and 6.3% of the optical and 
optical -i-NIR samples, respectively, are excluded from analyzes that require the stellar masses. 

We tested our stellar masses by comparing to those derived by a group at Max Plank us- 
ing SDSS photometry alone. Q The overall agreement is excellent. We find a modest offset in 
.SPIDER - M*MPi of 0.14 dex; this is because we use Kron magnitudes rather than fiber magni- 
tudes; Kron magnitudes are ~ 0.5 mag brighter, which translates into a difference ~ 0.2 in log{M). 
The rms scatter between the two mass estimators is only 0.1 dex, a factor of ~ 2 smaller than most 
systematic errors and biases introduced by different choices of IMF, extinction law, etc. More de- 
tails of the SED fitting and potential sources of errors and biases in the stellar mass estimates will 



'These mass estimates can be found at http://www.mpa-garching.mpg.de/SDSS/DR7/Data/stellarmass.htmll 



Table 1 . BC03 Galaxy SED Parameters 
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E(B - V) 
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be presented in an upcoming paper in this series devoted to this issue alone (Swindle et al. 2010, 
in preparation). 

We note that both the SED fitting and the analysis of color gradients (Sec. |7]) are performed 
with a difl'erent set of models (i.e. BC03) than that used for the spectral analysis (i.e. MILES, see 
Sec. 13.31) . The BC03 models provide a high resolution of 1 A in the wavelength range 3322 A to 
9300 A (griz), with a median resolution of 10 A in zYJHK. We do not use the MILES library here 
as we do for other tests in this paper due to its limited wavelength coverage (3525 A to 7500 A). 
With 7500 A lying near the central wavelength of the i-band, ignoring the izYJHK data inevitably 
yields failures in SED identification, provided the range of model parameters given in Tab. 1 . This 
uncertainty, coupled with the benefit of the K-band flux for mass-scaling, motivates this choice of 
models for the SED fitting. The technical justification for using the MILES library in analyzing 
spectra and color gradients is also detailed in Sec. 13.31 



3.2. Dynamical masses 

Dynamical masses, M^yn, are obtained from galaxy effective radii and velocity dispersions. 
Applying the virial theorem, we write 



M. 



dvn 



(1) 



where r^^K is the effective radius in K-band, and A is a scale factor that accounts for the internal 
structure of ETGs. Most previous works estimated A by assuming a single uni versal model for 
the light distribution in ETGs, such as the King model (|j0rgensen et al.lll996h. However, this 
does not account for the variety of ETG light profile shapes (e.g. ICaon. Capaccioli. & D'Onofrio 
1993b . as parametrized by the S ersic index. To account for this, following a similar procedure 
to that of iD'Onofrio et al.l (|2008i) . we estimate A using de-projected Sersic galaxy mode l s. Th e 
models and the de-projection procedure are detailed in ICoppola. La Barbera. Capacciolil (|2009h . 
The models describe spherical, non-rotating, isotropic galactic systems. We also assume that the 
spatial distribution of dark matter follows that of the stellar component. Under these assumptions, 
each model consists of a single Sersic component c haracterized by three parameters : the half- 
mass radius, the total mass, and the Sersic index (see Coppola. La Barbera. Capaccioli (|2009 ) for 
details). The use of more complex models (i.e. those with two-components) is well beyond the 
scope of this paper, and we postpone it to a forthcoming contribution. For each ETG, we construct 
the corresponding Sersic model, using the effective radius and Sersic index of the galaxy in K- 
band. We estimate A by computing the ratio of the quantities M^yn and cr^ x re^K for the given 
model. We compute ctq by projecting the model in two-dimensions, within a circular aperture of 
radius 1.5", i.e. the same aperture as the SDSS fibers. Since we do not have K-band structural 
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parameters for the optical sample of ETGs, Mdyn is only estimated for galaxies in the optical+NIR 
sample. Fig. [T] shows the dynamical masses versus the Petrosian magnitudes, from which we see 
that our limiting magnitude corresponds to a lower mass limit of ~ 3 x IO^^'Mq. 



3.3. Stellar population properties 

As described in Paper I , we have re-measured vel ocity dispersions for ETGs in the SPIDER 



sample using STARLIGHT dCid Fernandes et al.ll2005|) . This software determines the linear com 



bination of Simple Stellar Population (SSP) models which, broadened with a given cr, best matches 
the observed galaxy spectrum. Hereafter, we refer to the set of SSP models provided as input to 
STARLIGHT as the SSP basis, and to the output linear combination of SSPs as the synthetic 
spectrum. For the purpose of remeasuring the (To's, we run STARLIGHT with a basis of SSPs 
from the MILES galaxy spectral library (Vazdekis et al. 2010), covering a wide range of ages and 
metallicities, but with fixed solar abundance ratio ([a/Fe]=0). 

Here we also use STARLIGHT to estimate the Age, [Z/H], and [a/Fe] of each ETG. In con- 
trast to the more common practice of using index-index diagrams to estimate SSP-equivalent stellar 
population parameters |^ STARLIGHT uses all of the information present in the data, by fitting a 
set of SSP models to the entire galaxy spectrum. A detailed comparison of stellar population pa- 
rameter estimates from different techniques for the SPIDER sample is currently under way, and 
will be the subject of a forthcoming paper in this series. Here, we provide only the essential infor- 
mation on how we use STARLIGHT to derive Age, [Z/H], and [a/Fe]. For each galaxy, we run 
STARLIGHT using the alpha-enhanced MILES (hereafter or-MILES) spectral library (Cervantes 
et. al. 2007). As detailed in App.|Al the a-MILES library extends the MILES galaxy spectral li- 
brary to SSP models with non-solar abundance ratios, in the [a/Fe] range of -0.2 to -1-0.6, allowing 
the spectra of ETGs to be modeled with better accuracy than that achieved with solar-abundance 
MILES models. The a-MILES SEDs cover the same spectral range (3525 - 7500 A), with the 
same spectral resolution (2.3 A), as the (solar abundance) MILES library. Hence, they are well 
suited to analyzing SDSS spectra, whose spectral resolution is ~ 2.36 A FWHM, in the wave- 
length interval 4800-5350 A. To run STARLIGHT, we select a basis of 176 a-MILES SSPs, with 
ages ranging from 1 to ~ 14 Gyr in eleven steps, metallicities of [Z/H] = -0.68, -0.38, 0., 0.2, and 
enhancements of [a/Fe] = +0.0, +0.20, +0.40, +0.60. This grid in Age, [Z/H], and [a/Fe] allows 
us to cover a wide range of stellar population properties without overly increasing the execution 



^ We refer to SSP-equivalent parameters as those of the SSP model that best-matches the position of a given galaxy 
in the index-index diagrams. 
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Fig. 1. — Dynamical mass, for ETGs in the optical+NIR sample, as a function of the r-band Pet- 
rosian apparent magnitude, nip, from SDSS. Considering the scatter in the relation, and assuming 
the sample to be complete down to nip ~ 17.8, we conclude that the ETG's sample is complete, 
with respect to Mdyn, only above ~ 3 x lO^^M©. This mass limit is marked by the horizontal dashed 
line in the plot. The vertical and horizontal error bars in the lower-left of the Figure correspond 
(from left to right) to 25-, median, and 75-percentile errors on M^yn and nip. 
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time required to run STARLIGHT |j. For each galaxy, we first smooth the basis models to match 
the wavelength-dependent resolution of the SDSS spectrum (see Paper I for details). STARLIGHT 
outputs the light percentage, X, of each SSP basis model in the output synthetic spectrum. Hence, 
for a given property Y (e.g. Age), one can estimate its luminosity-weighted value, Y^, as 

Y - ^^^^ (2-, 
Yl - -Y^, (2) 

where the summation is performed over all the SSP basis models. The Age is obtained directly 
from Eq. |2l Since metallicity and enhancement are logarithmic quantities, we first apply Eq. |2] 
by setting Yi = 10^^^^^ and Y^ = 10^"^^^\ and then compute the logarithm of Y^ in both cases. 
The uncertainties on Age, [Z/H], and [a/Fe] are computed by comparing the estimates of these 
quantities for 2,313 galaxies with repeated observations in SDSS, following the same procedure 
described in Paper I to obtain the errors on ctq. The average uncertainties are ~ 20% on Age, and 
~ 0.05 for [Z/H] and [a/Fe]. 

An important potential bias when using the spectra from SDSS is the effect of a fixed spectro- 
scopic fiber aperture on the Age, [Z/H], and [a/Fe] estimates. The fraction of galaxy light inside 
the fiber will vary with galaxy size and redshift. Because we select the SPIDER sample based on 
the eClass parameter, the galaxies are selected to have passive spectra in their centers. This means 
that, as discussed in Paper I, our putative early-type sample could be contaminated by early spiral 
systems, i.e. spiral galaxies with a prominent bulge component. Indeed, in Paper I we show that 
this contamination is about 15%, falling to 5% for a higher quality sample defined on the basis of 
visual image classification. As we discuss at the end of Sec. 8.3 (Fig. 16), the trends with stellar 
population properties remain unchanged when restricting the SPIDER sample to the lower con- 
tamination subsample. The converse might also hold, i.e. we might be missing ETGs with strong 
emission features localized to their cores. However, very few ETGs are expected to have spectra 
exhibiting such strong central emission. 

In Fig. Owe display how Age, [Z/H], and [a/Fe] vary with the fraction of light inside the 
fiber wrt the total light, L(r< rfiber)A-tot, along with the distribution of this fraction, showing that 
~65% of the galaxies have a ratio between 0.2 and 0.6. The median trends, negligible compared to 
the scatter in Age, [Z/H], and [a/Fe], show that the derived stellar population parameters are not 
affected by the aperture effect, reinforcing what we have presented in App. |Bj i.e. that the average 
does not change much wrt the stellar population parameters, implying that the correlation of 
color gradients with stellar population properties has no dependence on rg. 



¥or 1,600 (out of 39,993 SPIDER) ETGs, we found that all SSP models in the basis, but those with [a/Fe] = 
+0.0, received zero- weight in the STARLIGHT synthetic spectra. For these objects, we enlarged the basis by including 
also a-MILES SSPs with [a/Fe] = -0.20, and re-ran STARLIGHT accordingly. 
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The values of Age, [Z/H], and [a/Fe] used here are not SSP-equivalent parameters, as com- 
monly used in the literature, but are luminosity-weighted quantities inferred from the galaxy spec- 
tra. As mentioned above, a detailed comparison of SSP-equivalent and luminosity- weighted pa- 
rameters will be performed in a forthcoming contribution. 



4. Color gradient estimates 



Internal color gradients are estimated using the same approach as La Barbera & de Carvalho 



(120091 hereafter LdC09), using the two-dimensional fits to the grizYJHK galaxy images with 
seeing-convolved Sersic models (Paper I). Color gradients are defined as the logarithmic slope of 
the galaxy radial color profile, Vvi'-x = dW - X)/d[logp), where X and W are any two wavebands, 
and p is the distance to the galaxy center. From grizYJHK data we can measure seven independent 
color gradients. As in LdC09, we calculate the color gradients between the g band and all other 
wavebands, in the form of Vg-x- For each galaxy and for each band, we estimate the mean surface 
brightness of the de-convolved Sersic model on a set of concentric ellipses, whose ellipticities and 
position angles are fixed to those from the r-band Sersic fit. The ellipses are equally spaced in 
equivalent radius p by 0.0 Irg ,. , where ,. is the the r-band eff"ective radius. The color index g - X 
at a given radius p is obtained by subtracting the mean surface brightnesses at the corresponding 
ellipse. Each colo r profile is fitted in the radial range of p^m = O.lr^. to p^ax = (see e.g. 



Peletier et al.lll990|) using an orthogonal least squares fitting procedure. The slope of the best- 
fitting profile gives the gradient ^ g-x- For the optical sample of ETGs, only the Vg_^, Vg_,, and 
^ g-z color gradients are estimated. As discussed in LdC09, small changes in the radial range used 
in this computation do not significantly change the color gradient estimates. 

To estimate the uncertainties on ^ g-x we repeated the above procedure by shifting and n in 
each band according to their measurement errors, taking into account the corresponding covariance 
term (Paper I). We repeat this using N = 1000 shifts, and the errors on V^.^ are then obtained from 
the widths of the distributions of Vg_x. 



5. Distribution of color gradients 



Fig. [3] plots the distributions of g - X color gradients for both the optical and optical -l-NIR 
samples of ETGs. For each distribution, we c ompute the corresponding peak and width {jx and 
cr) using bi-weight statistics (|Beers et al.lll990() . The values of yu and cr, with the corresponding 
bootstrap uncertainties, are reported in Fig. [31 as well as in Tabs. [2] and [3] for the optical -i-NIR 
and the optical samples of ETGs, respectively. All the distributions are sharply peaked. The peak 
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Fig. 2. — Correlation between Age (lower left), [Z/H] (upper left), and [a/Fe] (upper right) and the 
fraction of light inside the fiber wrt the total light, L(r< rfiber)/Ltot- The histogram of the fraction is 
exhibited in the lower right panel. The solid line displayed in each panel corresponds to the median 
behavior. 



-13- 



value of Vg-x smoothly decreases from about -0.07 in g - r to -0.3 in ^ - ^, i.e. the typical 
color gradient of ETGs becomes increasingly negative as we enlarge the wavelength baseline over 
which it is computed. It is remarkable that the peak values of Vg_,-, Vg_,, and V^.j., as well as the 
corresponding widths, are consistent between the optical (i.e. Tab.© and optical+NIR (i.e. Tab.O 
samples of ETGs. Fig. |3] also shows that the distribution of ^ - r through g - z color gradients is 
consistent between the optical (solid curve in the Figure) and the optical+NIR samples. The width 
of the Vg-x distributions smoothly increases from g - r through g - K. This trend might be due to 
larger errors on the near- infrared structural parameters than on the optical ones (Paper I), implying 
that the typical uncertainty on color gradients is larger for the optical-NIR rather than optical- 
optical Vj_x- Tabs.|2]and |3]provide the median uncertainties on the color gradients, cr'^'''', obtained 
by propagating the errors on structural parameters (Sec. H]). We estimate the intrinsic width of the 

CT' - ((T"''y] , with 

the constraint that cr' > 0. The corresponding intrinsic scatters are reported in Tabs. |2] and [31 It 
should be noted that this procedure is only a simplistic way of estimating the internal scatter, as it 
holds only for Gaussian distributions, which is not the case for the color gradient distributions (see 
Fig. 3). The cr' increases from g - r through g - K, implying that the distribution of optical-NIR 
color gradients has a significant intrinsic dispersion. The origin of this dispersion is analyzed in 
Sec. [81 where we study the dependence of stellar population gradients in ETGs on photometric 
and spectroscopic properties. We also note that the fraction of galaxies having inverted (positive) 
color gradients becomes less and less pronounced when we move from optical-optical to optical- 
NIR wavebands. For V^^^, we find that only a small fraction of ETGs (~ 8%) exhibit positive 
color gradients (V^.^^ > 0). From the peak value of ^g-x (-0.304) and the intrinsic width of the 
Vg_^ distribution (0.187, see Tab. 2), we estimate that only ~ 5% of galaxies in the optical+NIR 
sample are expected to have positive g - K gradients. Hence, the presence of a (small) fraction of 
objects with positive gradients is a real feature, and not merely explained by galaxies with negative 
gradients being scattered into the positive region by the measurement errors. 

The peak values of V^-x are consistent, within the uncertainties, with those estimated by La 
Barbera & de Carvalho (2009) (hereafter LdC09) for the entire SPIDER sample of 5, 080 ETGs 
with optical+NIR data available. A discrepancy is found for Vg„j, for which LdC09 measured a 
peak value of 0.089 + 0.006, while we measure -0.111 + 0.006. This difference does not affect 
the results of LdC09, since our value of -0.111 is fully consistent with those expected from the 
best- fitting stellar population models of color gradients in LdC09 (see their Fig. 1). The main 
differences between the sample of optical+NIR color gradients in this work and that of LdC09 are 
that (i) we select only the ETGs with better and more accurate structural parameters, and (ii) 
the computation of color gradients in LdC09 was slightly different than that performed here. We 
derive the color profiles using, for each band, the two-dimensional Sersic model whose axis ratio 
and position angle parameters are those measured in that band. In LdC09, we assumed the two 
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dimensional Sersic models in all bands to have the same axis ratios and position angles as those in 
the r-band. This also explains why the color gradient distributions are broader than those measured 
in LdC09 (see their Tab. 1). 



6. From color to stellar population gradients 

The color gradient reflects the radial variation of stellar populations in a galaxy. We can 
combine different estimates of color gradients from various pairs of available wavebands to infer 
the underlying stellar population gradient in an ETG. The approach we follow is presented, in 
general terms, in Sec. 16.11 la Sec. 16.21 we describe how this approach can be applied when the 
variation of stellar population properties is mainly due to age and metallicity. We discuss the role 
of the filters available to SPIDER in reducing the age-metallicity degeneracy. The results of this 
section are then applied to the color gradient estimates in Sec. |7]to constrain radial variations of 
age and metallicity within ETGs. 



6.1. Measuring variations of stellar population parameters 

We consider two stellar populations characterized by different physical parameters, such as 
age, metallicity, a-enhancement, etc. In general, the problem is how one can use the difference 
between the color indices of the two stellar populations to infer the corresponding difference(s) 
in physical parameters. The color index difference may be: (i) radial - the internal color gradi- 
ent; (ii) between galaxies with different masses - the Fundamental Plane coefficients in different 
wavebands (see La Barbera et al. 2010b; hereafter Paper II); (iii) in luminosity - the slope of the 
color-magnitude relation; or (iv) in the residuals in luminosity about a given scaling relation. We 
first describe the general approach we adopt, and then, in Sec. U\ we focus on its application to 

Table 2. Statistics of color gradients for the optical-l-NIR sample. 
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Vg-K 
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Fig. 3. — Distribution of optical+NIR color gradients of the SPIDER ETGs. From left to right and 
top to bottom, the histograms show the distributions of V^_x color gradients, with X = rizYJHK, 
for the optical+NIR sample of ETGs. On the y-axis, we report the fraction / of galaxies in the 
different V^^x bins. The peak and width of each distribution are reported in the upper-left corner 
of each panel, together with the corresponding Icr uncertainties. The vertical dashed lines mark 
the position of the peaks. The solid black curves in the three top panels show the distributions of 
Vg-r, Vg_/, and V^.^ color gradients, for the optical sample of ETGs. 



- 16- 



point (i). 



We indicate as A, the variation of a given physical property between the two stellar popula- 
tions, where the index i runs over the considered set of properties. The two stellar populations have 
differences in colors 6j, where j spans the set of colors defined by the available wavebands. One 
can introduce the following approximations: 

Si = y.A 



A,- 



(3) 



where A,y is the partial derivative of the j-th color index with respect to the i-th stellar population 
parameter. The above equation holds if (a) the color indices are well behaved functions of the stel- 
lar population parameters and either (bl) the A, are small enough to assume a linear approximation 
or (b2) the color indices are linear functions of the stellar population parameters. The quantities 
Aij have to be evaluated at the average values of the stellar population parameters of the two stellar 
populations. Eqs.[3]can be solved in ax^ sense, by minimizing the expression: 



(4) 



Requiring that the derivatives of with respect to Aj vanish, one obtains the following system of 
equations: 



(5) 



Setting Mjk = (Zr ^fa'^iv) ' (Zr ^y/^i) » we obtain a new system of N linear equations in the N 
unknown stellar population quantities Ay: 



z 



Mij-Aj=l 



(6) 



where the quantities Mij are completely defined by the partial derivatives A,y and the measured 
color differences 5,. If the matrix Mtj is non-degenerate, we can invert the linear system of Eqs. [6] 
and derive a unique solution for the unknown quantities A, (e.g. differences in age, metallicity, ...) 



from the measured color differences Sj. 



Table 3. Statistics of color gradients for the optical sample. 
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6.2. Constraining age and metallicity variations 

In this section, we investigate how the color indices in SPIDER can be used to constrain 
variations of stellar population properties in ETGs, specifically the effects of age and metallicity. 
Since we analyze color gradients in the form ^g-x (Sec. Sj), we consider here color indices of the 
form g - X, where X is one of rizYJHK. The ^-band includes the 4000 A break at the median 
redshift of the SPIDER samples, making the g - X colors most sensitive to the effects of stellar 
population variations. 

For the optical sample, each ETG is characterized by three independent colors, g - r, g - i, 
and g - z. Assuming that the relevant stellar population parameters are the age t and metallicity Z, 
equations [2] can be re- written as 



dig - r) ^ _^ d{g - r) 

d\og t ' d\og Z 

d{g-i) ^ d{g-i) 

d\og t d\og Z 



Sg-r = a,,, / •^+ r,,^ • (7) 



5g-i = ^TTT^ ■ + ^^rrrr • (8) 



d{g - z) , d{g-z) 
°s-z = ^1 — - • + — — • Az. (9) 
mog t dlog Z 

where A, and Az are the logarithmic differences of age and metallicity between two stellar pop- 
ulations. As noted above (see beginning of Sec. 16.11) . the quantities 6g-r, Sg-t, and dg-, could be 
the differences of color indices between galaxies having different parameters (e.g mass), although 
here we focus on the case where 6g^x is the difference in color index per radial decade in a galaxy, 
i.e. the internal color gradient (Sg-x = ^g-x)- In this case, the quantities A, and A, provide the 
difference in age and metallicity per radial decade (i.e. A, = V, and A-. = V^., see Sec. 17.11 ). Using 
the formalism of Sec. 16. 1[ A, and Az are derived by solving the linear system: 

M,A + M,zAz = 1 (10) 
Mz,At + Mz,zAz = 1 (11) 

where the matrix M is computed directly from the partial derivative of color indices with respect 
to log? and logZ, and the measured color differences 6g^r, dg-i, and 6g-^. For the optical -l-NIR 
sample of ETGs, the system of Eqs. |7]- [9] includes four extra equations, corresponding to the g-Y, 
g - J, g - H, and g - K color indices. The quantities A, and Az are still obtained by solving the 
system of Eqs. [TOl - fTTl where the matrix M is defined by using the color differences in all available 
wavebands. 

Since the NIR light is far less sensitive to metallicity (through line-blanketing) than the optical, 
and is also less sensitive to age, we expect that including the NIR wavebands helps to (partly) break 
the age-metallicity degeneracy, i.e. the degeneracy of the matrix M is reduced when including 
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all the g - X color indices. To examine this, we set 6i = 0.01 as a typical uncertainty in color 
differences. We adopt this value of 5, as it is an upper bound for the typical uncertainty on the 
mean value of color gradients. In fact, the typical error on the color gradient of a given galaxy 
is ~ 0.2-0.3 (Tab. 2). For the optical+NIR sample of ETGs, this implies a typical uncertainty on 
the mean color gradient of ~ 0.003-0.005 0. In order to estimate the derivatives of color indices 
with respect to age and metallicity, we adopt a polynomial approximation of the magnitude of a 
given stellar population model with respect to t and Z. Details of this approximat ion are provided 



in App. O Here, we consider SSP models from the iBruzual & Chariot! (|2003h synthesis code 
(hereafter BC03) with a Scalo IMF. We compute the matrix M for different values of t and Z of 
these SSP models, in the range where the polynomial approximation holds, i.e. 5 < t(Gyr) < 13.5 
and 0.2 < Z/Zq < 2.5. Fig. |4]plots the determinant of the matrix M, q = det{M), as a function of 
log t for different metallicities. The case where only optical colors are used is shown by the dashed 
curves, q is close to one for non-solar metallicity, while it is always smaller than one (by almost 
an order of magnitude) in the case of Z = Zq. This implies that for solar metallicity the matrix 
M is nearly degenerate and hence, as expected, optical colors alone are not able to effectively 
constrain age-metallicity variations. The solid curves show the cases where all the optical+NIR 
colors are used, q is then always greater than ~ 10. In particular, for solar metallicity, q is two 
orders of magnitude greater than with optical colors alone, implying that using all of the grizYJHK 
wavebands places much better constraints on age and metallicity variations. 

To illustrate this point more quantitatively, we perform Monte-Carlo simulations, setting A, = 
0.05 and Az = -0.3, which correspond to two stellar populations with a difference in age of 
~ 11% and a factor of two difference in metallicity. These values of A, and Az are close to 
those of the mean age and metallicity internal gradient of ETGs measured by LdC09. Using 
Eqs. |3l we compute the color differences, 6,, between the two stellar populations for all the seven 
g - X colors. The quantities A,^ are computed for t = lOGyr and Z = Zq. For each iteration, 
the color indices are shifted according to a normal distribution with a width of erg, and then we 
estimate A, and Az with the procedure outlined above. The procedure is repeated twice, adopting 
(i) OS = 0.002 and using only the griz colors, and (ii) with os = 0.005 and all the grizYJHK 
bands. The above cr^'s represent the typical uncertainties on the mean internal color gradients in 
the optical and optical-i-NIR SPIDER samples of ETGs ( see the errors on ft reported in Tab. |3] 
and Tab. [2l respectively). Fig. [5] shows the distributions of inferred A, and Az. In both cases, 
there is a strong correlation between A, and Az, implying that even with the full grizYJHK filter 



These values are obtained by dividing the typical errors on color gradients (~ 0.2-0.3), by the root square of 
the size of the optical+NIR sample of ETGs. Notice that the resulting expected errors on color gradients are slightly 
smaller than those on the peak values of ^g-x {~ 0.004-0.007, see Tab. 2), as the peak values are estimated by the 
bi-weight (rather than the mean) statistics. 
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set, one cannot completely break the age-metallicity degeneracy. However, the amplitude of the 
correlated variation is much smaller in the case where the optical and NIR bands are adopted. In 
fact, although for case (ii) erg is larger than in case (i), the corresponding error bars on A, and Az 
are much smaller, only ~ 0.01 for both t and Z. In other terms, optical+NIR colors constrain the 
size of the age and metallicity differences more effectively. From FigjH we conclude that even a 
small age gradient of ~ 11% can be detected at ~ 2.5cr with the optical+NIR dataset. 



7. Stellar population gradients 

In this section, we present our approach to constraining the variation of stellar population pa- 
rameters with distance from the ETG center. We consider only the effects of age and metallicity for 
simplicity, as the proper characterization of a given stellar population may require other quantities 
(e.g. the shape of the IMF). In App.lDl we describe the role of dust extinction gradients in ETGs, 
discussing their possible impact on our results. We first describe how we estimate the internal age 
and metallicity gradients (Sec. 17.11 ). In Sec. 17.21 we show how our approach allows us to define 
an effective color gradient, V^, that results from a combination of all the available color gradients. 
The systematics in V^^ estimates are then discussed in Sec. 17.31 



7.1. Inferring the age and metallicity gradients 

We adopt the procedure described in Sec. 16.21 to infer the radial gradients in age (V,) and 
metallicity (Vz). Setting 6g^x = ^g-x in Eqs.|7]-|9l the quantities A, and Az become the logarithmic 
differences in age and metallicity per radial decade inside a galaxy, i.e. A, = V, and Az = V^. 
Eqs.|7]-|9]are then rewritten as follows: 

- ^biT ^^ + ^1^ ^^' ^^^^ 

with X = rizYJHK. As shown in Sec. 16. 2[ in order to derive V, and Vz, one has to solve the linear 
system: 

M,,V, + M,zVz = 1 (13) 
Mz,,V, + Mz,zVz = 1 (14) 



where the matrix M is defined as in Sec. 16. 11 setting dj = ^g-x- The computation of M is performed 
using different stellar population models. For each model, we use Eqs.[T3]and [14] to derive the 
corresponding V, and Vz. This procedure, as presented in Sec. 6, minimizes the age-metallicity 
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Fig. 4. — Degeneracy, q, in deriving age and metallicity variations from color indices as a function 
of the logarithmic age of SSP models from the BC03 synthesis code. Different colors correspond 
to different metallicities Z, as shown in the upper-right corner of the plot. Dashed and solid curves 
correspond to the cases where optical and optical+NIR colors are used, respectively. 
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Fig. 5. — Measured A, and Az from Monte-Carlo simulations, where one shifts color indices 
according to the typical uncertainties on color gradient estimates, and recovers the true A, and Az 
via the approach described in Sec. l6.2[ Different colors correspond to cases where optical (black) 
and optical+NIR (red) color indices are considered. The error bars mark the 1 cr standard errors 
for the means of A, and Az recovered from the Monte-Carlo realizations. The horizontal and 
vertical dashed lines mark the true values of A, and Az. The dash-dotted lines are obtained from 
least-squares orthogonal fitting of Az vs. A,. 
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degeneracy depending on how many color terms are available and how they track the most con- 
spicuous spectral features (e.g. g - r probes the 4000 A break, yielding important information on 
the age of the underlying stellar population). 



We consider SSP models from three different sources: BC03, iMarastonI (|2005l) (M05), and 
Chariot and Bruzual (2009, in preparation; CBIO). These models are based on different synthesis 
techniques and have different IMFs. The M05 model uses the fuel consumption approach instead 
of the iso chrone synthesis of BC03 and CBIO. The CBIO code implements a new AGB phase 



treatment (iMarigo and Girardill2007h . The IMFs are: Scalo (BC03), Chabrier (M05), and Salpeter 
(CBIO). Moreover, we also use a composite stellar population model from BC03 having expo- 
nential star formation rate (SFR) with an e-folding time of t = 1 Gyr (hereafter BC03r=iGyr)- The 
models are convolved with the grizYJHK throughput curves, and the polynomial approximation 
described in App.|C]is used to calculate the partial derivatives entering in M. The derivatives are 
computed assuming an age oft = 9.27 Gyr, corresponding to a formation redshift of z = 2.5 at the 
median redshift of the SPIDER sample (< z >median~ 0.07), and solar metallicity (Z = Zq). For 
each stellar population model and each galaxy, we perform = 1000 iterations where in each of 
them we perturb the observed color gradients according to the corresponding uncertainties. The 
errors on V, and Vz are then the widths of their distributions resulting from all iterations. Note 
that in the above procedure all of the models are computed at the median redshift of our sample 
(z = 0.0725), neglecting the (small) redshift range of the sample. In App. we show that this 
approximation does not affect our results. 



7.2. The effective color gradient V* 

As discussed in Sec. 16.21 even using the entire grizYJHK set of wavebands we cannot break 
the age-metallicity degeneracy. This is further shown in Fig. |6l where we plot the age versus 
metallicity gradients obtained with the above procedure, for the optical-i-NIR sample, with the 
BC03 SSP models. The Vz and V, estimates are strongly correlated. As shown from the blue 
dashed line in Fig.|6l obtained from an orthogonal least-squares fit of Vz versus V,, the correlation 
is well described by a linear relation: 

Vzocff-V,. (15) 

The correlation is very similar to that between Az and A, in Fig. |5l This is shown by the red dotted 
line in Fig. [6l which is obtained by a linear fit of Vz versus Vj with slope fixed to that of the red 
dot-dashed line in Fig. [51 The slopes of the blue dashed and red dotted lines differ by only 4%, 
implying that the strong correlation of V^ and Vz is due mainly to correlated errors on differences in 
age and metallicity, such as the one between Az and A, shown in Fig. [51 In other terms, individual 
Vz and V^ estimates are affected by large, correlated statistical uncertainties because of the age- 
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Fig. 6. — Correlation of metallicity and age gradients as estimated for the optical+NIR sample of 
ETGs, with the BC03 SSP models. The correlation results from the age-metallicity degeneracy of 
color indices. The dashed line is obtained from a least- squares orthogonal fit of the data, with its 
slope, a, being used to define the effective color gradient V^^. The dotted red line is obtained by 
fixing the slope of the best-fit line to that of the red line in Fig. |5l The ellipses plot the typical 
measurement error contours that correspond to the median, 25- and 75-percentile uncertainties on 



Vz. 
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metallicity degeneracy (see the error ellipses in Fig. (6]). However, as shown in Sec. 16. 2[ using 
optical+NIR data allows unbiased estimates of the mean Vz and to be obtained, provided that 
the number of galaxies is large enough to make statistical uncertainties on the means small enough 
relative to the absolute values of the differences in age and metallicity ( i.e. and V^j. This allows 
us to bin galaxies according to different properties and perform a meaningful comparison of the 
mean V, and Vz among different bins (see Sec.|9l). 

In order to avoid the correlated uncertainties on individual and Vz estimates, we can define 
an effective color gradient: 

V* = Vz-a-V,. (16) 

From a geometrical viewpoint, V* measures the distance of a given point in the (Vz, Vf) plane 
along the Vz direction to the line defining the age-metallicity degeneracy. By definition, the cor- 
related variation of Vz and V„ which is described by Eq.[T5l does not change the V^,. Hence, the 
effective color gradient is not affected by the age-metallicity degeneracy. For galaxies having no 
age gradients, Vi, reduces to the galaxy metallicity gradient. Moreover, since a is negative (see 
below), Eq.[T6]can be rewritten as 

V, = Vz + \a\ ■ V,. (17) 

This shows that when Vz= 0, V^ is directly proportional to the age gradient. More generally. 
Eg. [T7] shows that V* behaves like a color gradient. A bluer outer stellar population relative to the 
inner one in a galaxy (because of either a metallicity or age gradient) yields a lower V^,. This is 
further shown in Fig. Ul where we plot V^, as a function of all the color gradients (see Sec. H]) for 
the optical+NIR sample of ETGs. On average, V^ monotonically increases as a function of each 
V^_x- In essence, can be seen as a color gradient that takes into account all of the colors. 

Using V* provides three important advantages over traditional color gradients. First, because 
V* combines all of the color gradient, its statistical uncertainty is significantly smaller than that on 
any color gradient using one pair of passbands. Second, different samples (bins) of galaxies can be 
compared using just one quantity (V^), rather than separately for each color gradient. Finally, as 
we prove below (Sec. 17.31) . the effective color gradient is a completely model-independent quantity 
(in the same way as the observed color gradients), which is not the case for Vz and V,. 

For each stellar population model, we derive a via an orthogonal least-square fit of Vz vs. 
Vf for all galaxies. Then we calculate V^^ for each galaxy, using Eq. \T6\ The value of a changes 
between different models, as well as for a given model when using a different set of wavebands 
(i.e. optical vs. optical -I- NIR). For the optical-l-NIR wavebands, we obtain a = -0.943 (BC03), 
-0.663 (CBIO), -1.341 (BC03r=iGyr), and -0.713 (M05). For the optical wavebands, we obtain 
a = -1.280, -1.186, -1.756, and -1.007, respectively. 
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Fig. 7. — Correlation of the effective color gradient, V^f, with the g - X color gradients, Vg-x- 
From left to right, and top to bottom, the X ranges from g through K. For each panel, the red curve 
is obtained by median binning the wrt the corresponding color gradient, while the ellipses 
plot the typical measurement error contours that correspond to the median, 25- and 75-percentile 
uncertainties on V*. 
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7.3. Systematics in V* 

We examine how sensitive V^^ is to varying the set of utilized wavebands. This is shown 
in Fig. [81 where we compare the V^, measured with optical+NIR and only optical bands. These 
two values of are linearly related, indicating that optical data alone can still provide a reliable 
estimate of the effective color gradient. This is consistent with the fact that V^^ is not affected by 
the age-metallicity degeneracy (influencing the optical more than the optical+NIR gradients) and 
is also proportional to each of the color gradients. On the other hand, the slope of the correlation 
between (optical) and (optical+NIR) is significantly greater than one at 2.00 + 0.03. This 
reflects the reduced statistical uncertainty on V^^ when incorporating the NIR information. This 
is further shown in Fig. |9] where we compare the histograms of V^, obtained using optical and 
optical+NIR data. The width of the distribution is significantly larger when using optical wrt 
optical+NIR data. Also, a small off"set in the peak of the two distributions is present. These 
results are quantified in Tab. |4] where we r eport the statistics , i.e. the peak (jj.) and width (cr) 
as estimated using the bi-weight estimator (jBeers et al. I99( h. for different samples and stellar 



population models. In the same Table, we also report the median uncertainty, cr^'"'', as well as 
the intrinsic width, cr'^, defined in the same way as the quantity cr' in Sec. |5l From these data 
we conclude that for the optical+NIR sample, the peaks of the V* distributions are consistent for 
different stellar population models, i.e. the estimated from optical+NIR is model-independent. 
In contrast, when using only optical data, a significant systematic different is present between the 
peak values of different models, though this difference is quite small (< 0.02) wrt the absolute 
peak value of V^. Furthermore, in agreement with Fig. 13 we see that cr is significantly larger when 
using optical data (cr ~ 0.4) as opposed to using optical+NIR data (cr ~ 0.25), i.e. V* is affected 
by significantly larger uncertainties (by a factor of ~ 1.6) when using only optical data. Finally, the 
peak values of the optical and optical+NIR are slightly different, in the sense that the effective 
color gradients are on average more negative by ~ 0.06 for the optical+NIR than for the optical 
samples. The differences in fi might be due to the different shapes of the optical and optical+NIR 
distributions, rather than some intrinsic difference. In fact, as seen in Fig.|9l the distribution of 
V* (optical) is clearly asymmetric. 

The robustness of V*, Vz and V, with respect to the adopted stellar population model is illus- 



Table 4. Statistics of the effective color gradients. 



model 


optical+NIR samp] 


e 






optical sample 










IT 


o-r 






tr 




< 


BC03 


-0.283 ± 0.007 


0.246 ± 0.005 


0.163 


0.184 


-0.219 ± 0.004 


0.407 ± 0.002 


0.377 


0.152 


M05 


-0.286 ± 0.007 


0.246 ± 0.005 


0.164 


0.182 


-0.209 ± 0.003 


0.388 ± 0.002 


0.355 


0.157 


CBIO 


-0.296 ± 0.007 


0.254 ± 0.005 


0.181 


0.178 


-0.235 ± 0.004 


0.436 ± 0.003 


0.406 


0.160 


BCCl-ir=iayr 


-0.291 ±0.008 


0.257 ± 0.005 


0.162 


0.200 


-0.225 ± 0.004 


0.417 ± 0.002 


0.390 


0.147 
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Fig. 8. — Comparison of the effective color gradient derived using either optical or optical+NIR 
wavebands. The dashed red line shows the orthogonal best-fitting relation between the two 
estimates. The slope of the relation is reported in the lower-right corner of the plot. The horizontal 
and vertical dotted lines mark the peaks of the distribution of V+ for the optical and optical+NIR 
data, respectively. The ellipses plot the typical measurement error contours that correspond to the 
median, 25- and 75-percentile uncertainties on the optical+NIR V^t. The ellipses are constructed 
by assuming, for simplicity, that the covariance term of the errors is proportional to that of the data 
themselves. 
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Fig. 9. — Comparison of the distribution of effective color gradients, obtained using either optical 
(black) or optical+NIR (red) data. The peaks of the two distributions are marked by the vertical 
dotted lines. Note how the distribution becomes narrower when including the NIR wavebands. 
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trated in Fig.dOl where we compare the V's obtained using the optical+NIR sample of ETGs for 
BC03 and M05. We see that the scatter in V^^ is significantly reduced in comparison to that for 
Vz or V,. As noted above, is model-independent, while this is obviously not the case for Vz 
and V, . This gives extra support to using as a useful and robust representation of the radial 
gradients of stellar population properties in galaxies. 



8. Dependence of the effective color gradient on galaxy parameters 

From Table 4 we see that the stellar population gradients exhibit significant intrinsic disper- 
sion, motivating an analysis of the correlations of versus other galaxy parameters. We study 
correlations of V* with several photometric quantities (Sec. 18.11 ). diff"erent estimates of galaxy 
mass (Sec. 18.21) . and spectroscopic parameters like central velocity dispersion and stellar popula- 
tion indicators like age, metallicity and [a/Fe] abundance ratio (Sec. 18.31) . Because of the large 
statistical uncertainty on individual measures of stellar population gradients, we bin the data and 
analyze trends with the mean . In each case, we verified that changing the number of galaxies 
in each bin does not aff"ect our results. The results presented in the following subsections refer to 
computed using the BC03 SSP models. Similar findings hold when considering other stellar 
population models described in Sec. U\ 



8.1. Effective gradients vs. photometric properties 

We bin V^^ with respect to the following photometric parameters: r- and K-band effective half- 
light radii, (Rej and r- and K-band total magnitudes (°*'^A/^ and °°^A//f), obtained through a 
PSF-convolved Sersic modeling of the galaxy image; logarithmic median Sersic index (< log n >); 
median axis ratio ( < b/a >); the 6(4 parameter, which characterizes the galaxy isophotal shape 
(boxiness - a4 < and diskyness - a4 > 0, see Paper I); the optical-NIR color index (r - K) and the 
median PSF fitting Ctpsp' Fig- [HI)- For each galaxy, median values are computed using data at 
the available wavebands, i.e. griz and grizYJHK. We adopted the r- and ^-band radii and mag- 
nitudes as representative of the sizes and luminosities of galaxies in the optical and NIR spectral 
regions, respectively. The r - K color index is estimated using 2DPH0T total galaxy magnitudes. 
Each bin includes the same number of galaxies, i.e. = 200 ETGs for the optical-l-NIR sample 
and N = 2000 ETGs for the optical sample. 

Fig. [TT] shows the relations between V^^ and the photometric quantities listed above. Panels 
(a) and (c) show that V* changes significantly with the r-band effective radius and r-band total ab- 
solute magnitude in the sense that larger and more optically luminous ETGs have stronger (more 
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Fig. 10. — Comparison of effective color (left), metallicity(middle) and age gradients (right) de- 
rived by using M05 and BC03 SSP models. Each panel plots the gradient difference in the sense 
of M05-BC03 as a function of BC03 gradients. The red dashed lines are the orthogonal best- fitting 
relations. The ellipses plot the typical measurement error contours that correspond to the median, 
25- and 75-percentile uncertainties on the x-axis variables. The covariance terms of the errors are 
derived as for Fig. [81 The left inset panel zooms into the center of the distribution to make the error 
ellipses visible. 
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Fig. 11. — Dependence of effective stellar population gradients on galaxy structural parameters. 
Each plot has been obtained by binning the with respect to a given quantity p. The following 
quantities are plotted: the r-band effective radius (panel a), K-band effective radius (panel b), r- 
band absolute magnitude (panel c), K-band absolute magnitude (panel d), logarithmic median of 
the grizYJHK Sersic indices (panel e), median of the grizYJHK axis ratios (panel f), the median 
of the a4 parameters in the gri passbands (panel g), r - K color index (panel h), median of the PSF 
modeling in the grizYJHK passbands (panel i). For each panel, the circles are the median V-^ 
in different bins of p. Error bars are 1 cr standard errors on the medians. The circles are connected 
by solid lines. The ETG sample with optical -i-NIR data is plotted in black. For panels (a), (c), (g), 
and (i), the binning is also done for the optical sample of ETGs (red color). 
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negative) effective color gradients. This result is somewhat expected as a more negative V^^ implies 
that bluer (either less metal rich or younger) stars in the galaxy are preferentially distributed to- 
wards the galaxy periphery. We note that the trend with ° ^^M^ is weaker than that with Re , - When 
we restrict ourselves to only optical data the trends with radius and magnitude are less pronounced 
(see below). Moreover, V^, has no dependence on radius or luminosity when we use AT-band ( Pan- 
els b and d). This finding is of particular interest, as the NIR light follows more closely the stellar 
mass distribution than the optical light. It implies that for bright ETGs the effective color gradient 
does not vary significantly along the stellar mass sequence. A mild correlation exists between 
and Sersic index, as shown in Panel (e). For logn>0.75 {n>5.6) ETGs with higher n tend to have 
stronger effective color gradients, likely because galaxies with higher n have higher luminosity dii e 



to the (optical) luminosity-Sersic index relation of ETGs (|Caon. Capaccioli. & D'Onofriolll993|) . 
Panel (f) shows no trend whatsoever of V* with axis ratio while with the boxy/disky parameter 
^4, shown in Panel (g), a weak correlation is present especially when optical-i-NIR data is used. 
In panel (h) we show that a strong correlation exists between V* and total galaxy colors. For 
r - K < 2.9, galaxies with bluer colors also tend to have stronger gradients. This is not a spurious 
result arising from the correlation of galaxy colors and magnitudes. In fact, the color- magnitude 
relation implies that brighter galaxies should have redder colors. This, together with the trend of V^^ 
with ° "^M^, would produce a trend opposite to that seen in Panel (h). Moreover, as shown in Pa- 
per I, no significant correlation is found between total magnitudes and total (rather than aperture) 
galaxy colors. Finally, Panel (i) probes possible systematics in color gradient estimates related to 
the PSF modeling. It shows that is not correlated with the;^f^ from fitting the galaxy image with 
a PSF convolved Sersic light distribution. Hence, our PSF models are accurate enough to provide, 
on average, unbiased V* estimates. 

In summary, we find that the effective color gradient in (bright) ETGs do not correlate with 
(NIR) luminosities and galaxy radii, while a correlation exists with galaxy colors. The former 
finding is further investigated in the next section, where we analyze the trend of with galaxy 
mass. The correlation of with colors is further studied in Sec. 18. 3[ where we bin with 
respect to the stellar population parameters of ETGs. 

The lack of correlation between and R^x (Panel b of Fig. [TT]) deserves further com- 
ments. In Paper II we show that the ratio of optical (^-band) to NIR (^-band) effective radii, 
Re,g/Re,K decreases as a function of Re,K, with the largest galaxies having Re,glRe,K ~ 1 (see Fig. 7 
of Paper II). This variation implies that the slope of the Kormendy relation for ETGs increases 
fro m ^ through K. A similar result, at optical wavebands {g - r), has been recently reported 



by iRoche. Bemardi. and Hydel (I2010t) . Using the ratio of optical radii as a proxy for the color 
gradient, the decrease of Re,g/Re,K with R^ k implies that larger galaxies should have shallower gra- 
dients (see also Tortora et al. 2010), in contrast with the lack of correlation in Panel b of Fig. [TT] 
However, color gradients are determined not only from the ratio of effective radii but also from that 
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of the Sersic indices. In fact, the variation of the surface brightness profile between two wavebands 
is determined from both Re and n. Fig.[T2lplots the logarithmic binned ratio of g to A'-band Sersic 
indices, log Hg/nx, as a function of log R^jc. We find that ng/n^ decreases as a function of Re,K, 
i.e. larger galaxies are more concentrated in the NIR than lower-7?^ systems. This compensates 
for the trend of Re,g/Re,K with Rg K, making V^, nearly constant with NIR radius. Thus there is no 
inconsistency between the lack of correlation in Panel b of Fig. [TT] and the correlation shown in 
Fig. 7 of Paper II. 

As mentioned above, we find that color gradients exhibit a complex behavior with galaxy 
luminosity. When estimated from optical+NIR color gradients, becomes more negative as the 
optical magnitude becomes brighter. On the contrary, for the optical sample of ETGs, i.e. when 
is estimated from V^_r, Vj_,, and V^_^, this steepening disappears, and does not vary with 
Mr (Panel c of Fig. [TTI) . The reason for the different behavior of V^^, when estimated using optical 
vs. optical+NIR data is illustrated in Fig. [131 where we plot the median color gradient, Vg^x, as 
a function of the optical magnitude for all the available wavebands {X = rizYJHK). The trends 
of Vj-x with ""^Mr are modeled with second-order polynomials (dashed curves in the Figure). 
The most remarkable feature is that the color gradient variation with magnitude depends on the 
waveband. The Vg_,- gradient exhibits a double-valued behavior, becoming flatter in both more 
and less luminous galaxies. For Vg-^, the double-valued behavior disappears and the V decreases 
monotonically as a function of °°^Mr. At redder wavebands, the trend reverses, and Vg^x becomes 
a monotonically increasing function of '^■^^Mr. As shown in Fig. Ul the effective color gradient 
is proportional to each of the individual color gradients from which it is estimated. When using 
optical data alone, mainly reflects the behavior of Vg^^, ^g-t^ and Vg_^. As a result, V^^ exhibits 
no strong variation with luminosity, similar to the double-valued behavior of V^-^- When 
is estimated from optical-l-NIR data, the optical-NIR gradients dominate the trend, with V* an 
increasing function of " ^^A/^. The main conclusion is that the behavior of color gradients with 
luminosity depends on the waveband where the luminosity is estimated (Panels c and d of Fig. [TTI). 
as well as the wavebands used to estimate the gradient itself. This complex dependence may 
explain, at least partly, the discrepant results fo und in the literature a bout the relation between color 
gradients and luminosity (mass). For instance. [Peletier et al.l (|1990|) first reported a surprising lack 
of correlation between color gradients and galaxy luminosity, noting that brightest ellipticals do not 
exhibit l ess steep gradien t s, as e xpected if they were the debris of repeated mergers of lower mass 
systems. Ide Propris et al.l(l2005h also reported no correlation between the size of the color gradients 



Tamura & Ohta 



and galaxy luminosity. 
steepe r (more negative) in brighter galaxies. The same resu 
(| 19941) for the bulges of spiral galaxies. On the other hand. 



bund that color gradients in cl uster ETGs become 



t was found by iBalcel 



La Barbera et all (120051) found that 
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color gradients do not depend on galaxy luminosity, and lChoi et al.l (l2007h found that optical color 
gradients are essentially constant over a wide range of galaxy luminosity, with a weak trend of 
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Fig. 12. — Logarithmic ratio of g to ^-band Sersic indices as a function of logRg in K-band. The 
soUd Une connects the data-points obtained by median binning the distribution of logng/riK with 
respect to logi^g /^, with each bin including the same number (N=200) of points. The trend implies 
larger galaxies to have a profile shape more concentrated in the NIR than in the optical (i.e. lower 
n^/riK). The dotted line marks zero Sersic index ratio. Error bars denote the Icr errors on the 
median in each bin. 
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becoming flatter at both fainter and brighter magnitudes. Recently, iRoche. Bernardi. and Hyde 
hereafter RBHIO), using the ratio of g- to r-band efi'ective radii from the SDSS as a proxy 
for the color gradient, also found that ETGs with intermediate luminosity have stronger gradients 
than those at low and high luminosities, in agreement with the double-valued behavior we observe 
for V„ 



8.2. Effective color gradient vs. galaxy mass 

Fig.dHplots the median V^^ as a function of galaxy stellar (M^) and dynamical (M^j.,,) mass. 
The binning is performed as in Fig. [TTl with each bin having the same number of galaxies. To 
quantify the dependence of V^^ on mass, we perform linear least-squares fits of V^^ with respect to 

and M^yn, with the slopes reported in the figure. We find no trend of the efi'ective color gradient 
with mass. In particular, the slope of V* vs. M^yn is consistent with zero at the 1.5 cr significance 
level, while the slope of V+ vs. is consistent with zero at ~ 0.5 cr. For M^yn < 2 x 10' 'Mq, the 

becomes more negative than at higher mass, but this effect might not be real, as in this mass 
regime, the SPIDER sample is incomplete due to the r-band magnitude selection (see Fig. [T]). We 
verified that the trend with stellar mass remains unchanged when using diff"erent stellar population 
models (e.g. CBIO) to perform the SED fitting (Sec. [3TB. 



8.3. Effective color gradient vs. velocity dispersion and stellar population properties 

Fig. [T5]plots the median efi'ective color gradient, , as a function of central velocity disper- 
sion, age, metallicity, and or-enhancement, for both the optical and optical -l-NIR samples. For the 
optical, the plot of V* vs. logcro seems to suggest a double- valued behavior, with V^, reaching 
a minimum around logcro ~ 2.2. However, this is not confirmed by the optical-i-NIR data, where 
no trend appears. This result supports the findings of Sees. I8.1l and l8.2[ that V* is not driven by 
galaxy mass. On the contrary, we find a strong dependence of V* on stellar population properties, 
i.e. the age, [Z/H], and [a/Fe]. In order to quantify these trends, we perform a robust linear fit of 
V* as a function of each of the three parameters, minimizing the sum of absolute residuals in V*. 
The uncertainties on fitting coefficients are estimated by the width of their distribution among 1000 
iterations, where the median V^, in each bin are shifted according to their error bars, and the fitting 
process is repeated. The slope and offset of the best-fit lines are reported in Tab. [5] From Panel (b), 
we see that the effective color gradient is flat for galaxies with older stellar populations. The trend 
is statistically significant for both the optical and optical-l-NIR samples, with the slopes of V^, vs. 
Age relations greater than zero at ~ 3 and ~ 9 cr, respectively. We note that galaxies with ages 
younger than ~ 4 Gyr seems not to follow the same trend as galaxies with older ages. As discussed 
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Fig. 13. — Median color gradient, Vg^x, with X = rizYJHK, as a function of r-band absolute 
magnitude, ''^^M^. Error bars denote 1 cr standard errors on the medians. The dashed lines are 
obtained by an second-order polynomial fitting of V^^^ versus ^■'^'^M^, with Vg~x as dependent 
variable. Different colors correspond to different wavebands, from g through K (top to bottom), as 
indicated by the labels on the left part of the plot. For better displaying the different curves, the 
Vg_x are re-normalized to have median values of -0.1, -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, from g 
through K. The actual medians are given by the peak values in Tab. |2l 
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Fig. 14. — Correlation of the median effective stellar population gradient, V*, with stellar (black) 
and dynamical (blue) mass, for the optical+NIR sample of ETGs. Error bars denote 1 cr uncertain- 
ties on the median value in each bin. Each bin includes the same number of galaxies (N = 200). 
Dashed lines mark the best-fitting relations, with the corresponding slopes being reported in the 
lower-left comer of the plot. The dotted lines define 1 <t confidence contours on the best-fit lines. 
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below, this inversion in the trend of V+ with Age is caused by galaxies with faint spiral-like mor- 
phological features and less accurate structural parameters. The behavior of V^^ with metallicity 
is less clear. For the optical sample, the slope of the V^, vs. [Z/H] best-fit relation is larger than 
zero at over 7 cr. However, the trend is only weakly detected in the optical-i-NIR sample, where 
the slope is only marginally larger than zero, at ~ 2.2 cr. We find that V+ strongly depends on 
the a-enhancement parameter, in the sense that galaxies with larger [a/Fe] also have shallower 
effective color gradients. The slope of the best-fitting line is significantly different from zero, at 
more than 5 cr, for both the optical and optical -l-NIR samples. 

As discussed in Paper I, the SPIDER sample of ETGs is affected by a small contamination 
from late-type galaxies with a prominent bulge component. To analyze the impact of this contam- 
ination, we define a purer sample of ETGs, where contamination is reduced to at most 5 %. To 
analyze if the trends of V^^ with stellar population properties are affected by the presence of late- 
type contaminants, we define two subsamples of ETGs, by matching the optical and optical-l-NIR 
samples with the lower contamination sample of Paper I. We also select only galaxies with better 
quality structural parameters, by removing objects whose two-dimensional fitting in any of the 
available wavebands is larger than 1.25 (see Paper I). The optical and optical-l-NIR subsamples in- 
clude 3, 928 and 31, 523 ETGs, respectively. Fig.[T6]plots the same trends as in Fig.[T5]but for the 
better quality subsamples. The coefficients of the best-fitting lines of V^, vs. age, metallicity, and 
enhancement, are reported in Tab. [6l The better quality subsamples exhibit trends fully consistent 
with those obtained for the full samples. In the case of V^^ vs. age, we do not find that galaxies 
with younger ages (< 4 Gyr) have flatter gradients, as found for the full samples. 

A further issue is that of the fixed aperture size (1.5" SDSS fiber radius) where the stellar 
population properties are measured. Because of radial gradients in stellar population properties, 
this fixed aperture size might produce a correlation of the estimated A^e, [Z/H], and [a/Fe], pa- 
rameters with galaxy radius. Since the V^'s correlate with the (optical) Re (see Panel a of Fig.fTT]). 
the aperture effect might also bias the correlation of the V's with Age, [Z/H], and [a/Fe]. How- 
ever, as shown in App.|Bl the variation of along ctq. Age, [Z/H], and [a/Fe], is negligible with 
respect to the full range of Re values (Panel a of Fig. [TTI) . implying that the trends in Fig. [15] are 
not affected at all by the aperture effect. 

Table 5. Coefficients of the linear fit of V^^ versus stellar population properties. 



Property 


optical+NIR sample 


optical sample 




offset 


slope 


offset 


slope 


log Age {Gyr) 


-0.529 ± 0.058 


0.227 ± 0.063 


-0.714 ±0.050 


0.498 ± 0.055 


[Zim 


-0.321 ±0.007 


0.225 ±0.100 


-0.263 ± 0.004 


0.430 ± 0.057 


[alFe} 


-0.239 ± 0.016 


-0.733 ±0.125 


-0.174 ±0.010 


-0.785 ± 0.086 
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Fig. 15. — The same as Fig. [TT]but for the spectroscopic parameters. From left to right and top to 
bottom, the median V^^ are plotted as a function of the following quantities: the central velocity 
dispersion (panel a), the luminosity-weighted age (panel b), the luminosity-weighted metallicity 
(panel c), and the luminosity weighted [a/Fe] abundance ratio (panel d). The optical-l-NIR and 
optical samples of ETGs are plotted with black and red colors, respectively, as shown in the lower- 
right corner of panel (a). 
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Fig. 16. — The same as Fig. [15] but for the optical+NIR and optical samples of ETGs with lower 
contamination from galaxies with faint spiral-like morphological features, and lower value of the 
two-dimensional fitting in r-band < 1.5). 



Table 6. The same as Tab. [5]for the samples with lower contamination and better quality of the 

structural parameters. 



Property 


optical+NIR sample 


optical sample 




offset 


slope 


offset 


slope 


log Age (G>t) 


-0.465 ± 0.062 


0.159 ±0.067 


-0.539 ± 0.046 


0.325 ± 0.050 


[Z/H] 


-0.319 ±0.008 


0.121 ±0.101 


-0.246 ± 0.005 


0.269 ± 0.059 


[a/Fe] 


-0.244 ± 0.017 


-0.626 ±0.130 


-0.174 ±0.012 


-0.621 ±0.094 
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9. Age versus metallicity gradients in ETGs 

In order to separate the contribution of age and metallicity to V*, we restrict the analysis to 
the optical+NIR sample of ETGs, as optical data alone are completely ineffective at disentangling 
the two components (see Sec. 16.21) . First, we analyze the statistics of the distributions of age and 
metallicity gradients (Sec. |97T1) . Sec. l9.2| deals with the trends of Vf andVz with respect to different 
proxies of galaxy mass, as well as stellar population parameters. In Sec. 19.31 the same trends are 
analyzed separately for low- and high-mass galaxies. In App. we also discuss some possible 
sources of bias in the trends of age and metallicity gradients, such as the role of internal reddening 
gradients in ETGs. 



9.1. Statistics of age and metallicity gradients 



In Tab. Ul we report the peak of the distributions of age (yu^ ) and metallicity (ju^^) gradi- 
ents for all the stellar population models descr i bed in Sec. U\ In agreement w ith previous works 
(e.g. IPeletier et al]ll99d : lTamura & Ohtall2003[ I2OO4I : La Barbera etaPboOBl) . we find that color 
gradients imply the presence of a negative metallicity gradient in ETGs (Vz < 0), with the outer 
stellar populations being less metal-rich than the inner ones. The peak value of Vz spans the range 
of -0.45 to -0.3, depending on the stellar population model. In agreement with our previous work 
(LdC09), we find that the peak V, is significantly greater than zero, i.e. ETGs have on average 
younger stars in the center than the outskirts. V, ranges from 0.05 to 0.2, being significantly larger 
than zero for all models. Hence, as discussed in LdC09, the existence of a small, but significantly 
positive age gradient is a robust (model-independent) result. From the values of Vz and Vf we see 
that, for all models, metallicity is the main driver of color gradients. 



Table 7. Peak values of the distributions of radial gradients in age (V,) and metallicity (Vz) of 

ETGs. 



model 




My, 


SC03 


-0.401 ± 0.021 


0.130 ±0.023 


M05 


-0.417 ±0.017 


0.198 ±0.025 


CBIO 


-0.309 ± 0.010 


0.052 ±0.012 


BCQ3r=lGyr 


-0.456 ± 0.023 


0.134 ±0.019 



log (Jq [km/s] Age (Gyr) [Z/H] 




Fig. 17. — Median values of the age (blue) and metallicity (red) gradients in ETGs as a function 
of (a) central velocity dispersion, (b-d) age, metallicity, and [a/Fe], parameters, and (e-f) stellar 
and dynamical masses. The plots refer to the optical+NIR sample of ETGs. Error bars denote 
1 (T standard errors on the medians. For each bin of a given quantity, the median effective color 
gradient, V*, is also plotted in black color. Age gradients are multiplied by the factor a, defined 
in Sec. Ul so that, for each bin, the effective color gradient can be approximately obtained by 
summing up the age and metallicity gradients in that bin (see Eq.fTTl). The vertical dashed lines in 
panel (f) mark the M^y^ used to define the low- (M^i < M^yn < ^^2) and high-mass (M^y^ > M^i) 
subsamples of ETGs (see the text). 
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Fig. 19. — The same as Fig. [I7]but plotting only galaxies in the high-mass bin {Mdyn > 8.5 x 
IQi^Mo). 
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9.2. Trends of age and metallicity gradients 

Fig.[l7]plots the median Vz (red circles) and V, (blue squares) in different bins of stellar pop- 
ulation parameters (age, metallicity, and enhancement), velocity dispersion, stellar and dynamical 
mass. For each bin of a given quantity, we also show the median V^^ in that bin (black triangles). 
Vjf, Vz and V, are estimated using BC03 SSP models. Age gradients are multiplied by the abso- 
lute value of the a parameter (see Sec. |7]), so that, for each bin, the sum of age and metallicity 
gradients gives, approximateljo V^^ in that bin (see Eq.[T6l). This allows us to analyze the relative 
contribution of age and metallicity to the trends of color gradients wrt the different quantities: 



9.2.1. Trends with mass proxies 

For massive galaxies, i.e. logcro>2.2 (ctq ~ 160 km s'^), M^>6 x 10"' Mq, and M^jy„>9 x 
10^°Mq, the lack of correlation between and either ctq, M^, or M^y„ (see Sees. 18.21 and 18.31 ) is 
the result of opposing variations of metallicity and age gradients wrt such quantities. The Vz de- 
creases while the V, increases as the mass proxy increases. The two trends cancel each other in the 
resulting color gradients. The origin of this opposite behavior is visible in Fig. [131 More luminous 
(massive) galaxies have flatter optical-optical color gradients (e.g. V^.^) and steeper optical-NIR 
gradients than lower mass systems. Because of the age-metallicity degeneracy, a flatter V^.^ may 
be due either to a positive age gradient in concert with a negative metallicity gradient, or an in- 
trinsically weaker metallicity gradient. The latter case is excluded by the optical-NIR gradients, 
as this would make the ^g-K flatten with mass, while we observe the opposite trend. Fig. [l7] 
also shows that low- (relative to high-) mass galaxies exhibit different trends. For logcr()<2.2 and 
M^<6 X IO'^Mq, no significant correlation of Vz and V, with mass proxy is seen, while for low- 
Mdyn objects (< 9 X IO'^Mq) a more complex behavior is observed. To illustrate this, we define 
two characteristic dynamical masses, Mdi = 2.5 x lO^^M© and M^j = 8.5 x IO'^Mq, marked by 
the vertical dashed lines in Panel (f) of Fig. [17] For M^„i<^d2 down to M^y,, ~ M^i, there is 
no significant correlation of Vz and V( with M^yn, while for Mdyn<Mai the metallicity gradient 
tends to decrease, while the age gradient tends to increase, as Mdy„ decreases. However, in this 
mass regime the SPIDER sample becomes incomplete because of the r-band magnitude selection 
(see Sec. [2]). We find that all the above trends are essentially unchanged regardless of the stellar 
population models used to estimate the age and metallicity gradients (see Sec. U}- 



^For each galaxy, V* is equal to the sum of age and metallicity gradients, according to Eq.[T6] On the other hand, 
the values of V*, V,, and Vz in a given bin do not satisfy exactly Eq.[T6] as they are obtained by computing the median 
value of the distributions of V*, V,, and Vz in that bin. 
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9.2.2. Trends with age 

For galaxies older than ~ 5 Gyr, we find that Vz decreases, while increases as a function 
of the Age parameter. The trend of V, is stronger than that of Vz, making the color gradient, V+, 
increase with Age (see Sec. 18.31) . For the M05 and BCO'ir=YGyr models, the behavior of and Vz 
with Age is the same as for BC03 SSP models. On the other hand, for C510, we find no significant 
variation of Vz with Age, while V, still increases with the Age. This confirms that the trend of 
color gradients with age, fox Age>5 Gyr, are mainly because of V,. For galaxies with Age < 5 Gyr, 
the trends of Vz and V, reverse. Lower ages imply more negative (positive) metallicity (age) 
gradients. Even in this case, the V, dominates the trend, and the resulting color gradient increases 
as the Age decreases. This result holds for all the different stellar population models. 



9.2.3. Trends with metallicity 

For galaxies more metal-rich than [Z///] ~ -0.12, the Vz tends to decrease, while the V^ 
tends to increase as metallicity increases. The two trends cancel each other producing no variation 
of the median color gradient, V*, as a function of the [Z///] parameter, consistent with the V*- 
[Z///] slope reported in Tab. [6] for the optical+NIR sample of ETGs. The trends of Vz and 
Vf seem to reverse at low metallicity, i.e. [Z///]< - 0.12 , similar to that observed for Age. In 
this regime, the variation of color gradient with [Z///] is dominated by Vz, which becomes more 
negative as [Z///] decreases. We found the behavior of V^ and Vz with Z/H to be independent of 
the adopted stellar population model. 



9.2.4. Trends with the enhancement 

In contrast to the trends with age and metallicity, we do not find any significant variation of 
Vf with [a/Fe]. Fig. [17] reveals that the strong correlation of V* with enhancement is because of 
the metallicity gradient decreasing with {a/Fe]. 



9.3. Low- and high-mass ETGs 

Since the trend of age and metallicity gradients wrt different mass proxies depends on the 
range of galaxy mass, we bin the sample of ETGs according to M^„n and analyze the trends of 
Vz and V, in each bin. We define two bins, including low- (Mji < M,/y„ < Mcti), and high- 
mass {Mdyn > Mdi) galaxies, where M^i and M42 are defined above. Although the sample of 
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ETGs becomes incomplete at Mdy„<3 x lO^^M©, we adopt a slightly lower mass limit of Mdy^ ~ 
2.5 X IO'^Mq, as this is approximately the Mdy„ at which the trend of V; and Vz with mass changes 
behavior (see Panel f of Fig. [T71) . 

Figs. [18] and [19] are the same as Fig. \TT\ but plot only galaxies in the low- and high-mass 
bins, respectively. In order to quantify the trends of V,, Vz, and V* wrt a given parameter, we 
fit each trend with a linear relation, minimizing the sum of absolute residuals of median gradients 
(see Sec. 18.11) . The slopes and their uncertainties are reported in Tabs. [8] and [9] for the low- and 
high-mass bins, respectively. From Figs.[T§]and[T3 as well as Tabs. [H] and [9] we see that: 



- Low-mass. The effective color gradient, V^^, does not show any significant correlation with either 

(To or Mdyn, while it tends to marginally decrease with stellar mass, as the V^^-M* slope in 
Tab. [8] is less than zero (-0.12 + 0.05) at the 2.2 cr confidence level. No significant trend 
of Vz and V^ is detected wrt any mass proxy (Mdy„, M^, and ctq). On the other hand, for 
stellar population parameters, we find that V^^ strongly correlates with Age and [a/Fe]. The 
effective color gradient increases for larger Age, and becomes more negative as the {ajFe] 
increases. The slopes of the V i,-Age and V i,-{alFe] linear fits differ from zero at more 
than 4 cr, proving that the trends are highly significant. The V^^ also tends to increase with 
metallicity, but the slope of the V^f-[Z///] fit is larger than zero at only 2.2 cr. As noted 
above, the trend of V^^ with Age is mainly due to V„ while that with {ajFe] is driven by Vz. 

- High-mass. For high-mass galaxies, the effective color gradient does not vary significantly with 

Age, metallicity nor with any mass proxy. In fact, the slopes of the correlations involving V^^ 
in Tab. [9] are all consistent with zero at less than 2 cr. For the [a/Fe], there is a significant 
trend of decreasing as the enhancement increases. The slope of the V^-{alFe] linear 
fit (-0.84 + 0.2) differs from zero at over 4 cr, and is consistent with that obtained for low- 
mass ETGs (-0.68 + 0.17, see Tab. [8]) and for the entire sample (-0.73 + 0.13, see Tab. [5]). 
For all quantities, we find anti-correlated variations of Vz and V,, which cancel each other 
in the trends of V^. For [a/Fe], the metallicity gradient dominates, resulting in the strong 
correlation of V^^ and [ajFe]. 



From Figs.[I71[T8]and[T9l we see that positive age gradients in ETGs are more associated with 
high rather than low mass galaxies. At high mass, the age gradient strongly increases as a function 
cro, Mdyn, and to a lesser extent with stellar mass, while these trends are not observed for low-mass 
galaxies. However, for both low- and high-mass, V^ becomes significantly positive for galaxies 
older than ~ 7 Gyr. This implies that age gradients in ETGs are mainly driven by mass, and only 
secondarily by Age. Panel (b) of Fig. [17] also shows that for Age < 5 Gyr, the trend of V, with 
Age reverses: the yo unger the age,the m ore positive the age gradient becomes. Consistent with the 
downsizing picture dCowie et al.lll996|) . galaxies at Age < 5 Gyr are low-mass systems, as shown 
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by the inversion in the Vf-Age trend being observed only for the low-mass bin (Panel b of Fig.fTSl). 



Stellar population parameters, such as Age, metallicity, and enhancement, are k nown to corre 



late with proxies of galaxy mass, such as veloc ity dispersion and stellar mass (see e.g. lThomas et al 



20051 : iNelan et alj 120051: iGallazzi et al.ll2006l) . Since for high-mass galaxies we find that internal 
gradients correlate with both mass and stellar population properties, a natural question is if these 
trends are just because of the correlation among stellar population properties and mass. In order to 
address this issue, we apply a correction procedure that removes the correlations of each V quantity 
with mass from all the trends in Figs. [TH] and [13 The procedure is illustrated in Fig. [201 where we 
show how the trends of the V's wrt to [Z/H] are corrected, for galaxies in the high-mass bin: 



i. First, we model the Vs-M^iy^ trends with third order polynomials in Mj^,„ (dashed curves in the 

top panel); 

ii. We compute the median M^m in each bin of \Z/H] (mid panel). As expected from the [Z/H]- 

mass correlation (e.g. [Thomas et al.ll2005|) . the mass increases as [Z/H] increases; 



iii. In each bin of [Z/H], using the corresponding median M^yn, we calculate the expected values 

of V*, Vz, and V; from the best-fit polynomials; 

iv. The expected V's in each [Z/H] bin are subtracted from the measured V in that bin; the sub- 

tracted values are re-normalized, by suitable additive terms, to the medians of the V^, Vz, 
and V, distributions (bottom panel in Fig. [201 ). Comparing top and bottom panels in Fig. [20[ 
illustrates how the correction procedure (partly) removes the correlations of V^^, Vz, and Vf 
with [Z/H]. 



We apply this correction scheme for both low- and high-mass ETGs, modeling the Vs-M^yn trends 
(point i. above) independently for each case. We obtain the following results: 

-High-mass. Fig. [2T[ exhibits the correlations of V^, Vz, and V,, with diff"erent quantities, for 
high-mass galaxies (M^y,, > Mai), after removing the V's-Mjj,„ trends. By construction, the 
procedure removes the correlations of the V's with M^yn (Panel f). However, it does not 
impact the trends involving stellar population properties. This is because such trends are 
due to the V - ctq correlations (rather than V's-M^y,,). This is shown in Fig. [221 where we 
repeat the entire procedure by removing the V's-ctq trends. The V's-ctq corrections remove 
all the trends, with the remarkable exception of the correlations among V^^ (Vz) and the 
enhancement. A linear fit to the V*-[a/Fe] trend in Fig. [22]gives a slope of -0.8 + 0.2, still 
fully consistent with that of the uncorrected trend (-0.84 + 0.19, see Tab.[9|). 
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Fig. 20. — Example of the correction procedure adopted to remove the correlation among the V's 
and dynamical mass from the trends of V^t, Vz, and V,, wrt galaxy parameters (FigsfT8]and[T9l). 
Top panel plots the V's as a function of M^yn for high-mass galaxies. Different symbols are the 
same as in panel (f) of Fig. \T9\ Dashed curves are third-order polynomial fits to the data. The 
mid panel plots the median value of M^yn in different bins of metallicity, [Z/H]. Metallicity bins 
as the same as for panel (b) of Fig. \T9\ The bottom panel shows the results of the correction 
procedure: for each metallicity bin, the polynomial fits are used to estimate the expected V's in 
that bin. The expected V's are subtracted from the data-points in the top panel. The resulting values 
are re-normalized, by suitable additive terms, to match the median values of the V^t, Vz, and V^ 
distributions. 
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Low-mass. In the mass range of M^i to Mdi, there is no variation of V^t with velocity dispersion 
or dynamical mass. Hence, the correlations between gradients and stellar population param- 
eters can not arise from age, metallicity, and enhancement varying with mass. However, at 
a given ctq, the stellar population parameters are correl ated with each othe r, with age being 
anti-correlated with both metallicity and enhancement (IGraves et al.ll2009l) . To account for 
this, we re-derive the trends between the V's and other parameters by removing the corre- 
lations with the Age. The procedure is the same as that adopted to remove the V's-Mjy« 
and V's-ctq correlations for high-mass galaxies (see above). Fig. [23] shows the same plots 
as Fig. [m but after the V's-Age trends are removed. We see that this correction does not 
affect significantly the correlations of the V's wrt metallicity and enhancement. In particular, 
after the correction is applied, a linear fit of wrt Z/H gives a slope of 0.29 + 0.12, fully 
consistent with that of 0.33 + 0.15 reported in Tab. [8l For V* versus [a/Fe], the slope is 
-0.59 + 0.15, still consistent with that in Tab. [H] (-0.84 + 0.19). It follows that, for low- 
mass systems, all the correlations of internal gradients with stellar population parameters are 
independent of each other, i.e. the gradients exhibit a genuine correlation with each single 
parameter (in particular age and enhancement). 



10. Summary 

In this paper, we examined a sample of 4,546 ETGs with grizYJHK data available. We con- 
sider an ETG to be a bulge-dominated galaxy with a passive spectrum in its central region (within 
the SDSS fiber). As such, the results obtained here do not apply exclusively to ellipticals, but to 
SO galaxies and galaxies with a red central bulge (although the effect of late-type contaminants is 
expected to be negligible, see Sees. 3.3 and 8.3). However, there is no a priori reason to expect that 
the formation mechanisms and evolutionary pathways are the same for all sub-classes of what we 

Table 8. Slopes of the linear fits of median V; (col. 2), V+ (col. 3), and Vz (col. 3), as a function 
of different galaxy parameters. Each line in the table corresponds to a different parameter, as 
reported in col. 1 . The slopes refer to the sample of ETGs in the low-mass bin. 





V, 


V, 


Vz 


log CTQ 


0.18 ±0.25 


0.050 ±0.104 


-0.16 ±0.24 


AgeiGyr) 


0.59 ±0.25 


0.406 ± 0.091 


-0.10 ±0.21 


[Z/H] 


-0.02 ± 0.30 


0.338 ±0.151 


0.33 ± 0.27 


[a/Fe] 


0.16 ±0.37 


-0.681 ±0.165 


-0.48 ± 0.38 




-0.01 ±0.17 


-0.111 ±0.055 


-0.05 ±0.15 


Mrf^„[10"Mo] 


0.19 ±0.16 


0.009 ± 0.068 


-0.13 ±0.16 
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Fig. 21 . — The same as Fig. [T9]but correcting the V in each panel for the correlations with dynam- 
ical mass. 
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Fig. 22. — The same as Fig. [19] but correcting the V in each panel for the correlations with ctq. 
Note how the correction removes all the correlations but those with the c-enhancement. 
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call ETGs. 

We measured structural parameters in all wavebands and stellar population indicators from 
the spectra. We describe a new approach to constrain age and metallicity gradients from the color 
information. Also, we define an effective color gradient, V*, that uses all the color terms provided 
by the photometric data and investigate its dependence on several parameters like age, metallicity, 
mass and [a/Fe]. Although the estimate of stellar population parameters (Sec. 13.31 ) and that of 
stellar masses and age/metallicity gradients rely on different stellar population models (BC03 vs. 
MILES), we find that all trends shown here remain unchanged when using a variety of population 
synthesis models to infer the gradients and estimate stellar masses. In the following we summarize 
our main findings: 

1) We present a new scheme to quantify stellar population gradients from color gradients. We in- 

troduce an effective color gradient, V^, which reflects all contributions from the optical+NIR 
colors. V* correlates well with all color gradients measured independently (as expected) and 
is model independent, while Vz and Vt (i.e. the metallicity and age gradients) are not. 

2) For photometric properties, we find that V^, has a mild correlation with the Sersic index; does 

not change at all with the axis ratio; and strongly correlates with total galaxy color - systems 
with bluer colors tend to have more negative V+. Using optical+NIR data we also find that 
larger and more luminous ETGs have more negative V*'s, although if we restrict ourselves 
to only ^-band data, shows no dependence on radius or luminosity. 

3) For galaxies more massive than M^y,, > 8.5 x 10'° Mq, the age gradient increases, while the 

metallicity gradient decreases as functions of mass, age, metallicity, and enhancement. The 
trends cancel each other in the color gradient, with the exception of those for the enhance- 
ment: the color gradient decreases as a function of [a/Fe]. All trends are a consequence of 
the correlations of the V's wrt both velocity dispersion and [a/Fe]. These parameters are 
the main drivers of the internal age, metallicity, and color gradients in massive ETGs. 

Table 9. The same as Tab. [8]but for galaxies in the high-mass bin. 





V, 


V* 


Vz 


logO-Q 


1.51 ±0.24 


-0.193 ±0.129 


-1.53 ±0.31 


Age{Gyr) 


1.54 ±0.33 


0.229 ±0.132 


-1.05 ±0.30 


[Z/m 


1.26 ±0.42 


0.006 ±0.165 


-0.98 ± 0.37 


[a/Fe] 


1.00 ±0.54 


-0.756 ±0.190 


-1.36 ±0.45 


MJ10"Mo] 


0.26 ±0.13 


-0.020 ± 0.052 


-0.26 ±0.13 


Mrf,„[10"Mo] 


0.56 ±0.15 


-0.082 ± 0.060 


-0.55 ±0.15 
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4) For less massive galaxies, 2.5 x 10^" Mq < M^yn < 8.5 x 10^" M©, no strong correlation of 

age, metallicity, and color gradients is detected wrt mass. However, color gradients strongly 
correlate with stellar population parameters, and these correlations are independent of each 
other. 

5) In both mass regimes, there is a strong anti-correlation between the color gradient and a- 

enhancement, that originates from the metallicity gradient decreasing with [a/Fe]. This 
could result from the star formation and metallicity enrichment being regulated by the inter- 
play between the input of energy from supemovae, and the temperature and pressure of the 
hot X-ray gas in ETGs. 

6) In both mass regimes, positive age gradients are associated with old galaxy ages (> 5 - 7 Gyr). 

For galaxies younger than ~ 5 Gyr, mostly at low-mass, the age gradient tends to be anti- 
correlated with the Age parameter, with more positive gradients at younger ages. 

We have studied the correlation of color gradients in ETGs with intrinsic galaxy properties, 
using an extensive sample, with excellent photometric and spectroscopic data over a broad wave- 
length baseline, and state-of-the-art analysis methods. However, even with such tools and data at 
our disposal, it remains difficult to gain substantive insight into the details of how galaxies actually 
form. The diagnostic value of color gradients will only become fully apparent when model predic- 
tions become sufficiently specific to realize the full potential of the high quality observational data 
and methods developed to estimate color gradients. 



A. The or-MILES models 

The or-MILES SSP models are constructed as described in Cervantes et. al. (2007), using both 
empirical and theoretical stellar libraries. The empirical libraries, like MILES (Sanchez-Blazquez 
et al. 2006), contain spectra of stars in the solar neighborhood, but mostly lack bulge-stars, with 
non-solar abundance ratios. Hence, they do not allow a detailed spectral modeling of metal-rich, 
alpha-enhanced systems, such as bright ETGs. The a-MILES models complement the empirical 
libraries, in the non-solar abundance regime, with the synthetic library of Coelho et al. (2005), 
consisting of high-resolution synthetic stellar spectra covering a wide range of stellar atmosphere's 
parameters. The Coelho et al. (2005) library covers both solar and alpha-enhanced mixtures in a 
wide baseline, from 3000 A to 1.4 //m, superseding previous versions by Barbuy et al. (2003), 
in the wavelength range 4600-5600 A, and Zwitter et al.(2004), in the range 7653-8747 A. The 

o 

resulting or-enhanced models consist of SED's covering the same spectral range (3525 - 7500 A), 
with the same spectral resolution (2.3 A), as the (solar abundance) MILES library. The or-MILES 
consists of 1, 170 SSPs, corresponding to: twenty-six ages from 1 to 18 Gyr; five [Z/H] from 
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-1.28 to +0.2, and nine [a/Fe] abundance ratios from -0.2 to +0.6. From this set of models, we 
run STARLIGHT by extracting a sub-set of 176 SSPs (see Sec. [33]). 

According to our tests, these preliminary alpha-enhanced models are clearly more efficient 
in fitting the ETG's spectra wrt to the solar-scale counterparts. For a subsample of 1,000 ETGs 
in the SPIDER sample, we have compared the quality of the STARLIGHT spectral fitting when 
using either solar-abundance or or-MILES SSPs. The fitting quality is measured, for each galaxy, 
by using: (1) the reduced of the fitting, x^/^A.ff, where A^^^^^ is the number of A values used 
in the fit minus the number of fitting parameters; and (2) the percentage mean deviation, < dev >, 
over all the fitted pixels between the input and synthetic spectrum. The < dev > is computed by 
taking the mean of \ObsA - Syn^l/ObsA, where Obs^ and Syn^ are the fluxes of the observed and 
synthetic spectrum at a given wavelength A, over the spectral interval 4320-6800 A. We found that 
the reduced decreases by 0.1 when using the alpha-enhanced models. The improvement in fit 
quality is illustrated in Fig. [24l Left panel shows, as an example, the synthetic spectra obtained 
for a given ETG, in a small spectral interval including the Mgb and FeSllO features, when using 
either solar (blue) or a-MILES (red) models. Although both models well describe the Fe, only 
the a-MILES synthetic spectrum (red) is able to reproduce also the Mg region. Right panel of 
Fig.|23]compares the <dev> parameters between solar and alpha-enhanced synthetic spectra. The 
a-MILES models provide a mean percentage deviation smaller than that obtained from the solar 
models by a factor of ~ 2. 

As a further test, we compared the solar MILES models with the or-MILES models having 
[a/Fe] = (i.e., we compared the models which use the empirical stellar spectra with those 
using synthetic ones, in the regime where both models should give the same answer). Fig. |25] 
plots the distribution of luminosity-weighted [a/Fe], obtained by running STARLIGHT with a 
basis of a-MILES SSPs on different solar-MILES SSP models (see Sec. 13.31 ). We considered a 
total of 84 solar models with a variety of ages (in 21 steps from 1 to 17.8 Gyr) and metallicities 
([Z/H] = -0.68, -0.38,0., 0.2). For each solar model, STARLIGHT was run with a basis of 
220 a-MILES SSP models, covering the same range in age and metallicity as for the analysis of 
SDSS spectra (Sec. [33]), and [a/Fe] = -0.2, 0, +0.2, +0.4, +0.6. The distribution of luminosity- 
weighted [a/Fe] in Fig.[25]is sharply peaked around zero, with a width less than ~ 0.01, proving 
the consistency of solar- and a models for the case of solar abundance ratio. 



B. Aperture corrections and color gradients. 

Since the V^'s correlate with the (optical) Re (see Panel a of Fig. [TT]), the fixed size of SDSS 
spectra might affect the correlation of the V's with Age, Z/H, and [a/Fe]. To address this issue, we 
show in Fig.[26]the median galaxy effective radius in different bins of ctq. Age, Z/H, and [a/Fe], 
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Fig. 24. — Performance of solar-abundance and a-MILES models to model the spectra of bright 
ETGs. (a) The spectrum of one ETG is shown, as example, in black color over a small spectral 
interval, including both Mgh and FeSllO features (gray bands in the plot). The Mgb index is well 
reproduced only by the alpha-enhanced models (red), as also seen in the residual spectrum (lower 
part of the plot), (b) Comparison of the mean percentage deviation, < dev >, between the solar- 
abundance and a-MILES synthetic spectra produced by STARLIGHT for a subsample of 1,000 
ETGs in the SPIDER (see the text). The percentage deviation decreases by a factor of ~ 2 when 
using the a-enhanced models. 
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Fig. 25. — Comparison of solar MILES models with or-MILES models having solar abundance 
ratio ([a/Fe~\ = 0). The Figure plots the distribution of luminosity-weighted [a/Fe], obtained by 
running STARLIGHT on SSP models from the solar MILES library, with a variety of ages and 
metallicities (see the text). The central and width values of the distribution are reported in the 
upper-left corner of the plot. Notice how the distribution is narrowly peaked around zero, proving 
the consistency of both sets of models for [a/Fe] = 0. 
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for the optical+NIR sample of ETGs. Red and black colors correspond to r- and K-band effective 
radii. As expected, some correlations exist between the Re and the spectroscopic parameters. In 
particular, the effective radius increases in galaxies with larger metallicity and a-enhancement. 
However, one can notice that the range of variation in amounts to only 0.1 - 0.2 dex, i.e. much 
smaller than the range of ~ 1.5 dex spanned by the optical and NIR effective radii (see Panel a 
of Fig. [TT]) . In the case of {ajFe] (Panel d of Fig. [26l) . the optical Re changes from ~ 0.55 dex at 
{a/Fe] ~ to ~ 0.7 dex at {ajFe] ~ 0.26. From Fig.[TTJ we see that the corresponding variation of 
is completely negligible, implying that the aperture effect is not affecting at all the correlations 
exhibited in Fig.[T5l 



C. Derivatives of SSP magnitudes with respect to age and metallicity 

We describe here how the colors available for the SPIDER samp le can be parametr i zed in 



terms of age and metallicity of stellar population models. We adopt the iBruzual & Chariot! (|2003|) 
(BC03) models. 

Fig. [27] exhibits the magnitudes of a simple stellar population model, normalized to a to- 
tal mass of one solar mass in stars, as a function of its age, t, for the grizYJHK wavebands of 
the SPIDER dataset. The griz and YJHK magnitudes are given in the AB and Vega systems, 
respectively. In order to describe the colors of ETGs, we consider relatively old ages, spanning 
the range of 5 to 13.5Gyrs. Different metallicities are also displayed in the range of 0.2 to 2.5 
times solar (Zq). Fluxes are computed with the BC03 code, convolving the SDSS and UKIDSS 
throughput curves with the model SEDs. In the optical passbands, the magnitude is essentially 
a linear function of log t. This is consistent with the fact that the flux of a simple stellar popu- 
lation is expected to chan ge as r", whe re a depends on the metallicity of the stellar population 



and the slope of the IMF (lTinsleylll973h . In the NIR passbands, the power-law approximation is 
less accurate, particularly at the lower age end. In order to calculate the color derivatives with 
respect to age and metallicity, for each waveband we fitted the corresponding magnitudes with a 
two-dimensional polynomial in log t and log Z. As shown in Fig. |27l a polynomial of degree eight 
reproduces model magnitudes very accurately. The rms of the fits are smaller than 0.01 mag in all 
the bands. Estimating fluxes using a composite stellar population with exponential star formation 
rate (SFR) and an e-folding time of 1 Gyr as a function of age and metallicity provides essentially 
similar results. Given the polynomial coefficients for each band, the computation of the color's 
derivatives is straightforward for any combination of t and Z. 
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Fig. 26. — Logarithmic median effective radius as a function of central velocity dispersion (Panel 
a), age (Panel b), metallicity (Panel c), and a-enhancement (Panel d). Black and red colors refer 
to K- and r-band effective radii. Error bars denote 1 cr standard error on the medians. Each trend 
is modeled by a linear fit, shown as a dashed line. The dotted lines mark the +lcr uncertainties on 
the fitted lines. The slopes of the fits and the corresponding 1 a uncertainties are reported in the 
left-upper comer of each panel. 
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Fig. 27. — Magnitudes of simple stellar population models from BC03 as a function of their for- 
mation epoch, t. The models are normalized to a total mass of 1 Mq (one solar mass) in stars. 
The griz and YJHK magnitudes are computed in the AB and Vega systems, i.e. the photometric 
systems of SDSS and UKIDSS, respectively (see Paper I). Different colors denote different metal- 
licities. The color become bluer from higher to lower metallicity. Four metallicities are plotted: 
Z = 0.2,0.4, 1.0, 2. SZq. Dashed lines show the result of fitting magnitudes as a function of age 
and metallicity with two-dimensional polynomials of degree eight. 
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D. Possible sources of bias in the Vs. 

The existence of radial gradients of absorption indices in ETG's spectra (e.g. Mgb and //g) 
implies that internal color gradients are likely driven by a variation of age and metallicity inside 
gal axies, consistent with our assumption in Eq. [121 However, as argued by ISilva & Wise 



and lWise & Silval (119961) . radial gradients in internal reddening might also contribute significantly 
to the observed color gradients. In such case, one would naturally expect a correlation to exist be- 
tween the gradient and the amount of dust, with stronger gradients in more dust-obscured galaxies. 
Although a detailed analysis of the dust distribution in ETGs is certainly beyond the scope of the 
present paper, we can evaluate if internal reddening might affect significantly our findings on the 
correlation among gradients and stellar population properties (Sec. |9l). As a rough estimate of the 
total amount of internal reddening, we use the color excess, E(B - V), obtained for each galaxy 
by the SED fitting procedure (Sec. [8]). Fig. |28] plots the median V* for the different E(B - V) 
(see Tab. 1). A clear correlation is detected, in the sense that steeper gradients are found at higher 
extinctions. The trend is mostly because of galaxies at high E{B - V) (> 0.3), with these objects 
amounting to less than 10% of the entire ETG's sample. We conclude that internal reddening might 
contribute significantly to color gradients only for a minor fraction of galaxies. Also, we find that 
selecting galaxies with E(B - V) < 0.2, for which the variation of color gradients with internal 
reddening is insignificant, does not change at all the trends shown in Fig. [T71 

Another possible source of concern is that color derivatives in Eq. [Tulare estimated by stellar 
population models reshifted at the median redshift, z = 0.0725, of the SPIDER sample, hence 
neglecting the effect of k and evolutionary corrections on color gradients. In order to see if this 
approximation impacts our results, we select only ETGs from the optical-l-NIR sample in a nar- 
rower redshift range of z = 0.0625 to 0.0825, leading to a subsample of 1, 984 galaxies. For this 
subsample, the peaks of the V^, V,, and Vz, distributions amount to about -0.31, 0.126, and -0.43, 
respectively, fully consistent with that obtained for the entire sample (see Tab.O. Also, the trends 
shown in Fig. [T7]turned out not to be affected by this selection. 
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UKIDSS project is defined in Lawrence et all (I2OO7I) . UKIDSS uses the UKIRT Wide Field Cam- 
era (WFCAM; Casali et al, 2007). The photometric system is described in Hewett et al (2006), and 
the calibration is described in Hodgkin et al. (2009). The pipeline processing and science archive 
are described in Irwin et al (2009, in prep) and Hambly et al (2008). UKIDSS data have been an- 
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Fig. 28. — Medians of the effective color gradient, V^^, as a function of the amount of reddening, 
E{B - V), in ETGs. The plot is obtained for the optical+NIR sample of ETGs. Error bars denote 
1 (T uncertainties on median values. The numbers of ETGs for different E{B - V) are reported 
above the corresponding error bars in the plot. Most of the galaxies (91%) have E(B - V) < 0.2, 
for which the variation of V* with color excess is smaller than ~ 0.02. 
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